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Chapter 1 


Introduction 


I believe in intuition and 
inspiration. 

Imagination is more important 
than knowledge. For knowledge is 
limited, whereas imagination 
embraces the entire world, 
stimulating progress, giving birth 
to evolution. It is, strictly 
speaking, a real factor in scientific 
research. 

Albert Einstein 


The physicist needs a facility in 
looking at problems from several 
points of view. 

Richard P. Feynman 


1.1 Aim of this thesis 

I believe that quantum correlations (also known as quantum entanglement ) ac¬ 
quired a special, almost mystical status in the collective conscious of physics 
community. My impression is that a lot of physicist tend to think of entan¬ 
glement as some elusive and incomprehensible property of quantum states. It 
exists, they believe, in addition to other well understood physical properties of 
quantum systems, such as coherence or symmetry, but it is always placed on a 
different level. Such a viewpoint is understandable for various reasons. For ex¬ 
ample, it is evident form the study of the literature dating before around 1995, 
that the concept of entanglement have not yet penetrated the vast majority 
of fields of physics (with the exception of topics revolving around famous EPR 
paradox and Bell inequalities, see Sec. E2). Nevertheless, it is clear that the 
advancements in all those fields where, and still are, unhindered by the lack 
of explicit considerations of entanglement. A reasonable explanation for this 
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state of affair is that in fact, the entanglement is not really an issue of its own. 
Rather, it is a complex, non-uniform in its nature, construct made out of a col¬ 
lective of intertwined “normal” physical properties. In other words, by focusing 
on those standard physical properties one can explain all occurring phenomena 
and never notice that formally the entanglement was involved. 

Does this mean that entanglement is an empty concept and should be dis¬ 
carded? I believe the answer is negative. One of the greatest struggles with 
understanding quantum mechanics is its remoteness from the everyday expe¬ 
rience of “classical beings” such as us. For me and many others, the question 
which aspects of the theory can be understood in terms of “semi-classical” mod¬ 
els, and which are inherently “quantum” has always been the most interesting 
one. Originally, the concept of entanglement has been conceived to help mak¬ 
ing this distinction. Presently, it evolved beyond the scope of this purpose, 
especially in the field of quantum information, but still I believe it is the best 
starting point we have for this kind of investigation. 

The aim of this thesis is to demystify the entanglement: to find out what it 
is in terms of physical properties of the system. I try to steer away from formal 
mathematical approach which could easily become detached from the physical 
intuition. Instead, I focus the discussion on building this intuition. The final 
result is the classical model constructed to represent quantum states, similarly to 
the celebrated Bohr’s model of hydrogen atom. The kinematic properties of the 
model are given by correlations present in the state. Thus, this approach might 
provide a new vantage point to examine the entanglement from the perspective 
of classical concepts we are all accustomed to. 

This thesis is not meant to simply summarize scientific results that I ob¬ 
tained during my studies; these can already be found in my publications. My 
understanding of the subject of non-classical correlations matured alongside the 
various other projects I was involved in. The preparation of the thesis was a 
great opportunity to collect my thoughts on the subject and organize them into 
a cohesive whole, which I could now share with the Reader. 


1.2 Non-classical correlations 

Although the foundations of quantum and classical physics are much different, it 
is often difficult to pinpoint which features of a particular system are intrinsically 
“quantum”. A good example of a “quantum” behavior is a wave-particle duality, 
which is a consequence of the ability of particles to exist in superpositions of 
quantum states. On one hand, the wave nature of massive particles is manifested 
in a Young double-slit experiment mm, which shows their ability to interfere, 
just like waves on a surface of a pond. On the other hand, the photoelectric 
effect [5j is a proof of a particle nature of the electromagnetic field, which is 
a wave, but also it consists of individual particles - photons. Even though a 
system which is a wave-particle hybrid is without precedence in classical world, 
still the wave or particle side of quantum phenomena on their own are perfectly 
conceivable in terms of classical physics. 

The most clear-cut distinction between “classical” and “quantum” can be 
made for systems composed of many particles when the properties of the en¬ 
semble are determined by the correlations between the constituents. Among the 
most important and evident types of quantum correlations are those between 
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identical particles. Classical physics allows for tracking every particle with per¬ 
fect accuracy without altering the dynamics of the system. Accordingly, even 
though the particles can be identical, i.e. characterized by the same set of fea¬ 
tures like mass or charge, they are always distinguishable by the virtue of their 
trajectories. The situation is dramatically different for identical quantum par¬ 
ticles. In quantum world the concept of trajectory simply does not exist due to 
Heisenberg’s uncertainty principle. Even if we were able to determine exactly 
the initial position of each particle it would be indeterminate in the following 
instant. Since, the particles cannot be tracked, they cannot be distinguished. 
Consequently, the physical properties of the system cannot change if any two 
identical particles are interchanged. Hence, the quantum state of identical parti¬ 
cles must be described by a wave function which is symmetric or antisymmetric 
with respect to permutation of particles. It is well known that this requirement 
has a profound consequences both for atomic-scale phenomena as well as for our 
everyday life as we know it. For example, Pauli’s exclusion principle, which is 
a direct consequence of identical fermions, like electrons and quarks, being de¬ 
scribed by antisymmetric wave functions, is an undergrid of the structure of the 
periodic table and overall large-scale stability of matter. Another example are 
bosons, which include photons, hclium-4 and Cooper pairs. They are described 
by symmetric wave functions, and as a consequence, tend to “bunch” together 
in the same quantum state. This bunching leads to phenomena like superfuidity, 
superconductivity and Bose-Einstein condensation of ultra-cold atoms. Overall, 
the concept of indistinguishability is unique to quantum mechanics and it has 
no counterpart in classical physics. 

Indistinguishability is not the only type of non-classical correlations. The 
famous gedankenexperiment proposed by Einstein, Podolski and Rosen (referred 
to collectively as EPR) revealed, then thought of as paradoxical, feature of 
quantum mechanics where a pair of particles in a particular quantum state would 
exhibit non-local properties [1]. Shortly after, Schrodinger made an attempt to 
extract the essence of non-classicality of EPR-type states, which led him to the 
concept of entanglement [5], which is the keynote of this thesis. The formal 
definition of entangled state used commonly nowadays reads: 

Definition 1. The state g of N parties is entangled if it is not separable, that 
is, it cannot be written as 



( 1 . 1 ) 


where gis a density matrix of party n and pi ^ 1 for all i with YhiPi = 1- 


Here “parties” may refer to subsystems composed of particles, sets of degrees 
of freedom (e.g. spin and position of an electron) or even parts of configuration 
space of a single particle - so called modes. The physical interpretation of 
entanglement follows from the definition of separable state. The state which is 
not entangled, i.e. can be written in the form of RHS of E3D, can be prepared 
“classically”. Each summand of RHS of (HUD is a product state, which means 
that it can be initialized by an independent measuring devices operating on 
each party individually. Suppose that these devices are equipped with a switch 
with settings i, set up in such a way that the device assigned to party n yields 
a state g^ as a result of the measurement. The set of devices can now be 
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supplemented with random number generator which would supply each of them 
with a choice of setting i with probability pi. Hence, the state prepared with the 
use of procedure described above is separable. The correlations between parties 
are determined only by the random number generator, which can be chosen 
to be a purely classical device. Thus, any state which could not be prepared 
by such a procedure has to contain correlations between parties that are not 
classical and is said to be entangled. 

Entanglement, by the virtue of definition, successfully formalize the notion 
of non-classical correlations. However, this discriminative definition is also the 
source of the limitations of the concept. We do not know what entanglement is, 
we only know what entanglement is not. It would be very naive to think that 
entanglement is “uniform” and there is no room for different types of quantum 
correlations that could be classified by their properties. The other problematic 
issue spawned by the general scope of the definition is an ambiguity of the 
concept of “party”. In most cases the physical situation explicitly defines the 
parties involved. Nevertheless, if the definition is invoked recklessly it may 
lead to nonsensical conclusions, like for example, equating a superposition with 
entanglement, which in turn can be destroyed or created by simple rotation of 
the reference frameQ. I shall address these two critical points in the following 
sections. 


1.3 Entanglement of modes 

According to the definition of entanglement, the system can be partitioned into 
parties in an arbitrary way, as long as the basis has been chosen so that it has 
a tensor product structure. Consequently, any superposition of states can be 
rewritten in such a way that it can be formally considered as an entangled state. 
To illustrate this statement, consider a two dimensional quantum system. I 
choose an orthonormal basis {|t), |4-)} and make a formal assignment mimicking 
the second quantization formalism 

|t)o|lt)|CU) and ||)o|0 t )|U) (1.2) 

By performing this mathematical trick I managed to turn an arbitrary super¬ 
position of the basis states into entanglement of two parties 

I*} = a It) + P II) +* a |l t ) |0j,) + p 10 t ) |U). (1.3) 

1 Consider two dimensional harmonic oscillator with Hamiltonian 

fr _ ( Px , rnuj 1 2 x 2 \ ( Pi muj 2 y 2 \ 

H + + + 

Written in a product basis of x and y directions, the state of single excitation in the direction 
x is separable 

1^) — &x |0os) |0j/) = |1®) |0y) 

Here di = yj ^ -pi ) is an annihilation operator in the direction i = x,y. However, 
the same state viewed in the reference frame rotated around z axis about ^ angle is entangled: 

IV>) = III) |0j,) = (cos^a x / -sin^aj,/) |0 X /) |0„/) = -h (|l x /) |0„/) - |0 X /) |l y /)) . 
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In this case, the entangled parties are “modes” defined as an orthonormal states 
that can be occupied by the particle composing the system. However, such an 
entanglement can be erased or created at will by a simple change of basis. For 
example, if I where to change the basis of our system to {|’F), |'F- L )} (|\E ,J ") is 
a state orthogonal to |\F)) and make a corresponding mode assignment 

|) o |l\if) |) and ||0ip) | ), (1.4) 

the state |\F) would be a separable state of modes \F and ’I'- 1 . These observations 
suggest that the concept of entanglement is merely a mathematical curiosity 
devoid of any physical significance. Indeed, this is the case if the entanglement 
of modes is considered in “vacuum”, without a context. Usually, the physical 
backdrop of a experimental setting favors particular choice of a modes describing 
the system, thus eliminating ambiguity in the definition of parties. A perfect 
example is provided by, the EPR experiment [3], which I will now discuss in 
more detail focusing on the issue of the choice of parties. 

Consider a system composed of two qubits, for example, photons which can 
be vertically and horizontally polarized. The system is initialized in a so called 
EPR state and then one of the photons is sent to detector located at far left 
and the other one to the detector located at far right. Thus, the synchronous 
measurements of polarizations of photons at the left and right site are performed 
on a state ket given by 

|EPR) = 7 f( l Vl\^)r + \^)l\X)r)- ( 1 - 5 ) 

Here the subscript L ( R ) indicate a state of a particle measured at the left (right) 
detector and |^) (|o)) designates vertical (horizontal) polarization of a detected 
photon. The state m is entangled without the shadow of doubt. However, it is 
not clear what exactly are the parties involved, although the answer might seem 
obvious at the first glance. It is tempting to say that the party described by ket 
| ip) L is the photon which went left, and \ip) R is a state of the second photon, 
the one which went right. This interpretation has to be dismissed immediately, 
because photons are indistinguishable and labeling them as “the one which went 
left/right” is meaningless. Such labels can only be attached to modes associated 
with states of definite values of physical quantities measured by left and right 
detectors, which as a classical objects, are distinguishable. Hence, Eq. OD is 
a proxy, or a shorthand of notation, for a more precise formula 

I EPR) = |0l,+>) \Qr,$) |l-R,-o-) + |0 l,$) |U,4+) |1r,$) |0i{,«.)) = 

= + a L,++ a kt) i VAC ^ ( L6 ) 

Here I used second quantization language: az,/ R ±/^ is a bosonic annihilation 
operator removing particle from vertical/horizontal mode located at left/right 
detector and |VAC) is vacuum state for which cq |VAC) = 0 for all possible 
modes i. The left and right modes are explicitly chosen by the measuring devices 
and it has been predicted that the entanglement of these particular modes is 
a necessary ingredient for observation of non-local properties of the quantum 
mechanics [B]. Thought experiment proposed by Einstein, Rosen and Podolsky 
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and the measurements on EPR state has been realized in laboratory □ BED and 
it has been demonstrated that the entanglement of modes can lead to observable 
phenomena unconceivable by classical physics. In fact, most of advances in the 
field of quantum information are based on utilizing the entanglement of modes 
chosen by a properly designed measurement schemes, including: teleportation, 
quantum cryptography, and quantum computing algorithms. 


1.4 Entanglement of particles 

The dependence on a choice of basis has proven to be the main difficulty with 
physical interpretation of mode entanglement. I argued that the ambiguity of 
assigning modes as parties can be lifted by a physical context and measure¬ 
ments associated with it. This dilemma cease to exist if particles are designated 
as parties instead of modes. The Hilbert space of many particle system is a 
tensor product of subspaces describing each particle. Separable state of parties- 
particles cannot be made entangled, and vice versa , by means of the change of 
basis. Indeed, the choice of basis states cannot compromise or alter the identity 
of particles, hence the only allowed basis transformations are local to particle 
subspaces. Consequently, the question whether particles are entangled or not is 
independent of the choice of basis - a property very appealing form the physical 
point of view. 

The entanglement of distinguishable particles is no different then the en¬ 
tanglement of modes. As long as particles are not identical a separable states 
can be prepared classically, because in principle each party can be measured 
and initialized independently. Such a classical procedure of state preparation 
is impossible to implement for indistinguishable particles, since the measuring 
devices that could address parties individually cannot exist. Thus, the notion 
of equating separable states to classical states have to be reconsidered when 
dealing with identical particles. For starters, it might seem that because of 
symmetrization/anti-symmetrization condition enforced on a state ket all states 
of identical particles must be non-separable. In fact, this is the case for fermions 
but is not for bosons, to which I now turn our discussion. 

I begin by stating that a separable pure state of N identical bosons | (j>\N) 
must be a product of N identical single-particle orbitals | <j>) !10], i.e. 

\^N) = \^f N . (1.7) 

When the bosonic field operator acts on the state ( 11 . 71 ) . the result is 

\4>\ N) = VN(j>{x) \<j>\ N — 1), (1.8) 

which is a fixed -N counterpart of the property of a coherent state of light |<f>) 
defined by the relation £^ + \x) |$) = |$), where £^ + \x) is the positive- 

frequency part of the electromagnetic field £{x). In Eq. ( 11 . 81 ) . (j>{x) is a single¬ 
particle function determining the spatial properties of the system. Note that an 
analogical state of fermions is impossible and thus, particle-separable state of 
fermions cannot exist. In general, the separable state of N identical bosons is 
a mixture of states m , thus follows the definition of particle entanglement of 
bosons: 
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Definition 2. The state g of N identical bosons is particle entangled if it cannot 
be written as 



(1.9) 


Here T>(j) denotes the integration over complex field <f> and is a probability 
distribution, i.e. it is normalized and its integral with every F((f >) ^ 0 over any 
volume V is non-negative, 



( 1 . 10 ) 


v 


There is a direct analogy between the so-called P-representation of the state 
of light, where the density matrix is represented as pi; g ht = / 2?$P($) |<£} ($|. 
If functional P is not a probability distribution (i.e does not satisfy condition 
similar to (11.101) 1. the electromagnetic field is considered to be non-classical 
uni nu mi- Analogically for N indistinguishable bosons, if condition (11.101) is 
not fulfilled, the density matrix cannot be written as a statistical mixture of 
separable coherent state, meaning that correlations between the particles are 
genuinely quantum. 

On the other hand, the separable states of bosons are classical because the 
symmetrization enforced by the indistinguishability on state EB is only su¬ 
perficial. In principle, identical state can be prepared out of independent, dis¬ 
tinguishable particles by initializing each one of them in the same quantum 
state using a proper set of measuring devices. Hence, the indistinguishability of 
bosons in separable state is inconsequential and such a state can be simulated 
by a classical system; at least for as long as the processes which address parties 
individually are absent. 

An entanglement of identical particles is often considered as “useless” as 
opposed to mode entanglement which has proven to be a valuable resource 
for many applications in the field of quantum information. Many protocols 
utilizing entanglement are designed under the assumption that each party can 
be address individually, thus excluding the possibility of using indistinguishable 
particles. Nevertheless, as I already argued previously, correlations of identical 
particles are very important and have profound consequences. I shall discuss the 
usefulness of particle entanglement, which does not have to be considered only 
in terms of applicability of certain class of protocols, in the upcoming sections. 

1.5 Useful entanglement 

The difficulty in grasping the nature of entanglement is a direct consequence 
of a general and discriminative definitions © and ©• Possible solutions to 
this problem is to examine the subject form the utilitarian point of view. The 
idea is that certain tasks, like cryptography or computing, can be performed 
better when entangled states are used instead of states that are only classically 
correlated. This line of reasoning leads to the concept of classification of en¬ 
tanglement by the degree of their usefulness for a given task. Now I proceed to 
formalize this abstraction. 

Let (g) be the efficiency at which a given task T is performed using 
state g. Since the form of separable state is known it is often possible to find a 
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maximal efficiency achievable with non-entangled states: 

Q t (q separable) T Q(| 


( 1 . 11 ) 


Therefor, if a state allows for efficiency greater then the classical bound Qq~, 
then the state has to be entangled: 


Q r (g) > Q% => g is entangled. 


( 1 . 12 ) 


The reverse implication is not guaranteed, i.e. not all entangled states are able to 
outperform classically correlated states. A criterion such as (11.121) gives means 
to classify types of entanglement, in this case, in terms of their usefulness for 
a given task. Moreover, the efficiency Q can be treated as an indicator of the 
degree of non-classical correlations - the greater the Q the more of the useful 
entanglement is present. 

This approach to entanglement detection is reminiscent of the method of 
entanglement witnesses mm and criteria such as spin squeezing nananzmsi 
and Cauchy-Schwarz criterion m- The efficiency criterion ( 11 . 121 ) posses an 
advantage over other tests that, by its definition, it automatically provides an 
application for useful entanglement and the physical context. Therefore, the 
properties of usefully entangled states can be juxtaposed with everything that 
is know about the task it is useful for, which can turn out to be a valuable 
source of physical intuition. On the other hand, the entanglement detected by 
the other mentioned criteria can be characterized only by the fact that it is 
detected by those criteria. For these reasons I will adopt the notion of useful 
entanglement and the efficiency criterion as a basis for further discussion on 
subject of nature of non-classical correlations. 






Chapter 2 

Atomic interferometer 


In previous chapter I argued for the advantages of using the efficiency criterion 
(11.121) as a tool for investigating the properties of entanglement. However, in 
order to implement such approach, first one has to decide which task T is to be 
performed with the help of non-classically correlated states. My task of choice 
is quantum metrology , in particular the quantum interferometry. 


2.1 Quantum metrology 

Metrology is the science of measurement. Its main goal is to develop methods, 
both practical and theoretical, for measuring properties of physical systems as 
precisely as possible. Precise measurements are essential in our everyday life, 
without it we would not be able to build safe and efficient buildings, machines, 
medications and so on. Metrology is also of great importance for science itself. 
In the end any scientific inquire boils down to the measurement and more precise 
results allow for drawing more decisive conclusions which can help in testing a 
theory or point us in new directions of research. 

Any measurement procedure consists of three steps: the preparation of a 
probe, its interaction with the system to be measured, and the probe readout. 
However, this process is inevitably affected by statistical or systematic errors. 
The source of the former can be accidental (e.g. deriving from an insufficient 
control of the probes or of the measured system) or fundamental (e.g. deriving 
from the Heisenberg uncertainty relations). Whatever the origin of the error, 
it can be reduced by repeating the measurement and averaging the resulting 
outcomes. According to the central-limit theorem, given a large number N 
of independent measurement results each having a standard deviation a, the 
average error converges (with growing N ) to a Gaussian distribution with the 
standard deviation equal to a/y/N, so that the resultant error scales as 1 /\/N. 
On the other hand, a single measurement realized with a probe composed of 
N uncorrelated sub-probes (e.g. particles forming a matter-wave) includes the 
independent repetitions automatically thus yielding the error which also scales 
as 1/y/N. This behavior is referred to as the shot noise limit (SNL) and I 
will demonstrate in the upcoming chapter that it is also associated with proce¬ 
dures which do not fully exploit the quantum nature of the probe. This implies 
that SNL can be surpassed only when one employs quantum effects, such as 
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the entanglement among the probing particles utilized for the measurements. 
Consequently, the SNL is not a fundamental quantum mechanical bound as it 
can be overcome by using non-classical strategies. Quantum metrology studies 
the fundamental bounds on precision imposed by quantum mechanics and the 
strategies which allow for attaining them. More generally it deals with measure¬ 
ment and discrimination procedures that receive an enhancement in precision 
through the use of quantum effects. Therefore, the precision of measurement can 
serve as a efficiency criterion sensitive to “quantum effects” and in particular it 
can be used for detecting quantum entanglement between particles constituting 
the probe. 

Atomic interferometers |19j form a family of devices which exploit the wave 
nature of matter, providing a unique and powerful tool for modern quantum 
metrology. They have been used to measure atomic properties [2U], to study 
quantum degenerate systems ED, for precision measurements E2, and are very 
sensitive probes for inertial effects with applications in gravimeters and gyro¬ 
scopes EaUHEHUMj. 

Bose-Einstein condensates (BECs) are promising candidates for atom inter¬ 
ferometry owing to their macroscopic coherence properties. Following the first 
observation of the BEC interference in 1997, various building blocks of the BEC 
interferometers have been realized individually. The interference experiments 
with BECs were performed using Bragg beams in ballistic expansion, with freely 
propagating matter-waves in a guide EH Eg |22| . By splitting a single trapped 
BEC into two separated clouds in a double well, interference was observed after 
switching off the trapping potential EHlEDlSa- These and many other examples 
show that a long-standing goal of realizing a full interferometer with ultra-cold 
atom systems is in our reach. 

A fundamental difference between photon and matter-wave optics is the pres¬ 
ence of atom-atom interactions which might lead to generation of non-classical 
correlations. The unprecedented degree of control over these interactions via 
Feshbach resonance techniques [33. allow for preparation of a BEC probe in 
strongly entangled states which can beat the shot noise sensitivity of the inter¬ 
ferometric measurements. This is the main advantage of atomic interferometry 
over its optical counter-part: with BECs it is relatively easy to entangle thou¬ 
sands of particles, while with light the current technical limitations only allow 
for few photon entanglement. 

In the following section I describe a few important examples of experimental 
realizations of atomic interferometers based on BECs and introduce a theoretical 
tools for describing such systems. 


2.2 Examples of atomic interferometers 

The basic design of an interferometer involves two “arms” through which the 
probe can travel and the ability to “mix” the signals from the two ports in a 
coherent manner. In case of atomic interferometers the “arms” can be realized 
by well separated modes; this includes modes in the real space, momentum space 
or even in the space of internal degrees of freedom such as Zeeman levels. The 
beam-splitter operation which transfers the particles between modes as well as 
phase difference imprint can be achieved by properly devised external fields. 
Having this picture in mind I proceed to describe in more detail some of the 
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Figure 2.1: Example of atomic interferometer: double well setup m a) Ini¬ 
tially Bose-Einstein condensate is created in the optical lattice with atoms oc¬ 
cupying 2-6 well separated sites. The harmonic potential with tunable frequency 
is imposed onto the system. The potential is adiabatically “squeezed” forcing 
atoms in two adjacent sites, which define the arms of the interferometer, b) 
The height of the barrier between potential wells can be manipulated in order 
to control the rate of coherent tunneling, which allows for realization of beam¬ 
splitter operation. The relative depth of the wells can also be manipulated by 
introducing the energy offset and thus realizing the phase imprint operation. 


most important experimental realizations of atomic interferometers. 

The first paradigm is the double well setup achieved with condensate trapped 
in the optical lattice. In [53] the particles in BEC of s 'Rb atoms where dis¬ 
tributed over a small number of lattice sites (between two and six) in a one¬ 
dimensional optical lattice. The occupation number per site ranges from 100 to 
1,100 atoms. The two modes representing arms of the interferometer where the 
two states of the external atomic motion corresponding to the condensate mean- 
field wavefunctions in two, well separated adjacent lattice sites (see Fig. 12.11 (all. 
The beam-splitter operation which allows for mixing particles occupying the two 
modes was realized by coherent tunneling through potential barrier separating 
the lattice sites. By adjusting the height of the barrier it is possible to control 
the rate of tunneling or even stop it all together. The imprint of the phase differ¬ 
ence between arms comes from the difference in depth of lattice sites (see Fig. 12.11 
(b)). The difference could be caused by the variation of external potential (e.g. 
gravitation, magnetic field) on the length scale of wells separation thus allowing 
for high-resolution measurement of this variation. The information about im¬ 
printed phase can be inferred from absorption imaging measurements of atom 
number in each lattice site. The wells of the lattice can be fully resolved thanks 
to imaging techniques with a resolution of 1pm, which is well below the lattice 
spacing of 5.7pm. This allows for the determination of the atom number in each 
lattice site by direct integration of the measured atomic density. Local interfer¬ 
ence measurements after a condensate expansion time short enough that only 
neighboring sites overlap reveal the phase between these wells, confirming that 
the coherence of the condensate is preserved. The non-classical correlations be¬ 
tween particles where provided by the on-site repulsive atom-atom interactions 

m- 

A full interferometric sequence composed of beam-splitter followed by phase 
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imprint followed by another beam-splitter with trapped 8 ' Rb BEC confined on 
an atom chip was demonstrated in [35] . The interferometric scheme relied on the 
coherent splitting and recombination of a BEC in a tunable magnetic double-well 
potential, where the matter wave is confined at all times. Thanks to a spatial 
separation of ~ 2pm between the two wave packets, the geometry was sensitive 
to accelerations and rotations. By tilting the double well out of the horizontal 
plane for a variable time, the energy difference was applied and thereby imprint 
a controlled relative phase between the interferometer arms. A non-adiabatic 
recombiner translates the relative phase into an atom number difference, which 
is directly read out using a highly sensitive time-of-flight fluorescence detector. 
As in the previous case also here the particle interactions in BEC matter waves 
lead to a nonlinearity which generated the entanglement between atoms. 



Figure 2.2: Example of atomic interferometer: internal degrees of freedom m 

a) The Bose-Einstein condensate of atoms in F = 2 hyperfine state is confined in 
optical dipole trap at a homogeneous magnetic field. The spin dynamics in atom- 
atom collisions lead to creation of pairs in spin up and down state (tuf = ±1/ 

b) In order to read out the populations of the Zeeman modes, the magnetic trap 
is turned off and the strong magnetic field gradient is turned on. In the fashion 
of the Stern-Gerlach device, the internal modes of the atomic cloud become 
spatially separated and ready for density measurement with standard imagining 
techniques. 

Now I turn to a different paradigm where the arms of the interferometer are 
realized with modes of internal degrees of freedom. The experiments reported 
in [36] started by creating a condensate of 2.8 x 10 4 87 Rb atoms in the hyperfine 
state F = 2 with horizontal spin orientation (Zeeman substate tuf = 0) con¬ 
fined in an optical dipole trap. The spin dynamics in particle collisions where 
used to create up to 10 4 paired neutral atoms in spins up and down states 
(mj? = ±1). These collisions are bosonically enhanced if the output modes 
are occupied. Therefore, they act as a parametric amplifier for a finite initial 
population in tuf = ±1 or for pure vacuum fluctuations. During the paramet¬ 
ric amplification of vacuum, the total number of atoms produced in mj? = ±1 
and its fluctuations increase exponentially with time (see Fig. 12.21 (a)). The 
conjugate variable of the total number is the sum of the two atomic phases, 
whose fluctuations are exponentially damped. Furthermore, the number dif¬ 
ference between tuf = ±1 atoms is zero (without fluctuations), and hence the 
corresponding conjugate variable, the relative phase, is fully undetermined. The 
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underlying physics closely resembles that of optical parametric down-conversion 
in nonlinear crystals, currently the most important technique to generate non- 
classical states of light. The spin dynamics where initiated at a magnetic field, 
where an excited spatial mode is populated and vacuum fluctuations are am¬ 
plified. The states F = 2, ijif = ±1 are populated by spin dynamics for an 
optimal duration of 15ms. The internal-state beam-splitter was implemented 
by driving the transition connecting the F = 2, rriF = ±1 states with three res¬ 
onant microwave pulses. Subsequently, the dipole trap was switched off and all 
three spin components where recorded by absorption imaging after they where 
spatially separated by a strong magnetic field gradient (see Fig. 12.21 Ib'l'l. 



Figure 2.3: Example of atomic interferometer: twin matter beams a) 

Special geometry of the trapping potential creates the effective two-level system, 
b) An optimally controlled shaking of the trap inverts the population from the 
ground state to excited state with almost 100% efficiency, c) The excited state 
is metastable. Due to collisions pairs of atoms get demoted to ground state and 
the excess of energy is converted into motion with opposite momenta directed 
along x-axis due to appropriately shaped trapping potential. 


As a last example I refer to an experiment described in mi where it was 
demonstrate how collisional deexcitation of a one-dimensional degenerate Bose 
gas can be used to efficiently create matter wave beams separated in momentum 
space. The starting point is a quasi-BEC of 87 Rb atoms magnetically trapped 
in a tight waveguide potential with a shallow axial harmonic confinement (say, 
along a;-axis ) on an atom chip. The scheme relied on an effective two-level sys¬ 
tem in the radial vibrational eigenstates of the waveguide (see Fig. l2.3l faB. This 
was accomplished by creating unequal level spacings in the radial y, z-plane by 
radio frequency dressing, which introduces anharmonicity and anisotropy. Due 
to the increasing level spacings, the ground state and the first excited state 
along y, have the lowest energy difference among all possible combinations, es¬ 
tablishing a closed two-level system. Having prepared the gas in the ground 
state, the population was inverted by transferring the atoms almost entirely 
to the excited state (the efficiency of coherent transfer reached ~ 97%). The 
transition is driven by shaking the trap along the radial y-direction on the scale 
of the ground state size (~ lOOnm). The trajectory of the shake (total duration 
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5ms) has been optimized by an iterative optimal control algorithm (see Fig. 12.31 
(b)). In the experiment, the displacement was achieved by driving a current in 
an auxiliary chip wire, parallel to the main trapping wire. The population inver¬ 
sion represents a highly non-equilibrium state of the system, analogous to a laser 
gain medium after a pump pulse. For the ensuing relaxation, due to wave-guide 
geometry the only allowed channel was a two-particle collisional process, emit¬ 
ting atom pairs with opposite momenta. Within a binary collision, two atoms 
are scattered from the excited state to the ground state and the excess potential 
energy is transformed into the kinetic energy of each atom. Due to momentum 
conservation each one acquires momentum of the same magnitude but opposite 
direction along the elongated axis of the trapping potential (see Fig. 12.31 (c)). 
Similarly to the previous example the emission process can be understood as 
a matter wave analogue of degenerate optical parametric amplifier, where the 
initially empty twin-modes are seeded by vacuum fluctuations and gain an ex¬ 
ponentially growing population if the phase matching conditions are fulfilled. 
Finally, the beam-splitter can be achieved by applying a Bragg pulse which can 
coherently bring matter-wave into superposition of counter-propagating beams. 
The populations of twin-beams was measured with fluorescence images once the 
trap potential is switched off and the atoms propagate freely separating from 
the source. 

In the last two examples the source of twin-mode matter waves was a colli¬ 
sional phenomena taking place among atoms forming the condensate. Also, due 
to indistiguishability of atoms and bosonic enhancement, the scattering into well 
separated modes is responsible for formation of strong non-classical correlation 
among the particles. Reference [3B reviews the Bogoliubov theory in the con¬ 
text of such twin-beam experiments and describes the process of pair generation 
leading to highly entangled states useful for interferometry. The nature of this 
entanglement will be discussed in the upcoming chapters. 


2.3 Pseudo-spin of two-mode atomic interferom¬ 
eter 


According to examples presented in the previous section Bose-Einstein conden¬ 
sate undergoing interferometric measurement can be considered as an ensemble 
of N qubits defined by two distinct modes L and R constituting the arms of 
interferometer. Due to bosonic nature of atoms the state of the system for fixed 
N can be described in terms of occupation numbers of each mode. The mode 
occupation number states form an orthonormal basis in this N + 1 dimensional 
space and are defined with the help of annihilation and creation operators which 
add and subtract a particle in a given mode: 


Basis of 

mode occupation 
number states 



ri)R = 


(a[r ( 4 ) (JV - ra) 

Vn\ y/(N — n)\ 


| VAC) 


( 2 . 1 ) 

TO=0,l,...,iV 


Here a,L /r are the annihilation operators for mode L/R which satisfy standard 
bosonic commutation relations {{o-l/RjOcr/r} = 1 and all other commutators 
equal zero), |VAC) is a vacuum defined as a state which satisfies a_L/_R|VAC) = 0. 
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Incidentally, mode occupation number states are also eigenstates of operator 


f _ - a R a R _ til - flR 

z ~ 2 ~ 2 


This operator, together with 


f _ a^aR + a^ R a L 

X - 2 

d^dij - a ] R a L 
Jy ~ 2i 


( 2 . 2 ) 


(2.3) 

(2.4) 

(2.5) 


constitutes a set of three orthogonal components of fictitious spin operator. 
Straightforward calculation shows that this set of operators indeed satisfy com¬ 
mutation relations of angular momentum algebra [, J 3 ] = ieijkJk- The action 
of these operators on basis states is given by 


J z \n) L \N - n) R = ( nL 2 nR ) \n) L \N - n) R = [N 2 n) \n) L \N 


n) R 


+ jy “I" 

(' flL+flR ) 2 , 

4 ' 


n) L \N - n) R = 

n) L \N - n) R = y(t + 1 )l n )il^ r 


flL+flR 

2 



( 2 . 6 ) 

Therefore, the system composed of N bosonic qubits is equivalent to single 
pseudo-particle with spin N/2: 


{\ti}l\N ^) fl} n=0j ]JV ° { IAf )} j- n /2,m=-j,...,j 


(2.7) 


This formal equivalence can always be established as long as the two-mode 
approximation holds. Moreover, pseudo-spin operators can also be ascribed a 
transparent physical interpretation. The component J z is an observable asso¬ 
ciated with population imbalance between the modes of interferometer. As it 
was noted in the previous section, the high resolution imaging techniques allow 
for very accurate evaluation of the number of atoms in each mode and hence 
the population difference is the main property of the system accessible through 
direct measurement. The x and y component combine into ladder operator 
J + = J x + iJ y = a\a R and J_ = J x — iJ y = cl^cll which describe a process 
where a particle is transfered from one mode to the other. These operators are 
associated with a beam-splitter operation, which is a necessary ingredient of 
operational atomic interferometer (see previous section). 

The equivalence between the two-mode state of bosonic qubits and a pseudo¬ 
particle with spin will prove to be very productive in the upcoming chapter 
where I use it to analyze types of particle entanglement found in ultra-cold 
atom systems. 
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Chapter 3 

Efficiency of interferometer 


In this chapter I provide a basic introduction to the theory of parameter es¬ 
timation, on which the efficiency criterion for atomic interferometers is based 
on. 


3.1 Distinguishability of quantum states and the 
problem of parameter estimation 

The principle of operation of the atomic interferometer is based on detecting 
changes in a state of the probe induced by the interaction with the measured 
system. The evolution of the probe state obviously depends on the properties 
of the system; in particular it depends on the unknown value of the parameter 
9 characterizing its certain aspect. The goal is to estimate the exact value of 
this parameter. The efficiency at which this task is performed is tied to the 
precision of the estimation procedure. 

Once the interaction with the system is concluded we end up with the output 
state which has some information about the parameter imprinted onto it, 

interaction with 

A the measured system ^ / sw /„ _ N 

Pin ---> Pout (6>) (3.1) 

To gain access to this information a measurement (in a sense of quantum me¬ 
chanics) has to be performed on p 0 ut- Preferably, the measurement should be 
chosen so that it addresses those properties of the output state which where 
most affected by the parameter imprint. In the final step of the procedure the 
data acquired from the measurement is used to estimate the real value of 9. 

An alternative point of view is that the value of 9 parametrizes a path 
through the space of quantum states followed be the probe under the influence 
of the system (see Fig. 13.11) . Therefore, estimation of the value of 9 is directly 
related to the ability to pinpoint the position of the state on this path. Then, 
the precision of the estimation relays on the capability for distinguishing the 
neighboring quantum states occupying the path. The precision can, thus, be 
increased by preparing the probe in the state which is most susceptible to the 
evolution driven by the interaction with the system. Indeed, if the state un¬ 
dergoes a significant change when exposed to the system it should be easier to 
determine a small variation of the parameter. In other words, the state which 
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moves along the path with greater “speed” will travel a greater distance even 
for small increments of 9 (see Fig. 13.21) . 




Figure 3.1: The evolution of probe state under the influence of the measured 
system viewed as a 9-parametrized path in the space of density matrices. The 
ability to locate an output state g 0u t = q( 9) on the path is required for estimating 
the value of 8. 

This point of view was adopted in [39] where the authors formulated the 
problem of parameter estimation in terms of distinguishing neighboring quan¬ 
tum states in a space of density matrices. They have found measurements which 
optimally resolve neighboring states, and characterized their degree of distin- 
guishability in terms of a Riemannian metric, increasing distance corresponding 
to more reliable distinguishability. These considerations allowed to establish a 
type of uncertainty principles which allowed to relate the precision of estimation 
with the susceptibility of the probe state to change under process of parameter 
imprint. Below I present a brief summary of the derivation of this result. 

3.2 Distinguishability metric 

I begin by reviewing a derivation of the distinguishability metric for probability 
distributions [40] . After drawing v samples from a probability distribution, one 
can estimate the probabilities p(£) as the observed frequencies /(£), The proba¬ 
bility for the frequencies is given by a multinomial distribution, which for large 
v is proportional to a Gaussian a exp [ — (i//2 )(/(£) — p(£)) /p(0] ■ A nearby 
distribution p(£) can be reliably distinguished from p(£) if the Gaussian exp [ — 


(i//2)(p(£) — p(0) 2 /p( 0] is small. Thus the quadratic form (p(£) — p(f)) 2 /p(£,) 


provides a natural distinguishability metric on the space of probability distri¬ 
butions (PD), called the statistical distance: 



(3.2) 
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oifi+dQ) 




Figure 3.2: States which are more susceptible to the influence of the measured 
system can undergo a significant change and cover grater distances in the density 
matrix space even for very small variations of parameter 9. For such states it is 
easier to distinguish g{9) from the neighboring g{9 + d9) which allows to locate 
the state on the path with greater reliability. This translates directly to precision 
of parameter estimation. 

The notion of the statistical distance has been generalized by the authors to 
mixed quantum states and thus obtain a natural Riemannian geometry on the 
space of density operators. 

Consider now a curve g{9) on the space of density matrices. Performing a 
measurement on the state is the only way available in quantum mechanics to 
distinguish g{9) from the neighboring matrices g{6 + dQ). In general quantum 
mechanical measurements are described by a set of non-negative, Hermitian 
operators 11^ which are complete in the sense that 



(3.3) 


Such a set is refered to as positive operator valued measure (POVM). The quan¬ 
tity £ labels the “results” of the measurement; although written here as a single 
continuous real variable, it could be discrete or multivariate. The probability 
density for result £, given the parameter 6 , is 


pm = Tr(fl^(0)) . 


(3.4) 


Thus, the density matrix g(9 ) has been mapped onto probability distribution 
which can be treated with classical statistical distance dsp D : 



F(mA)do 2 , 


(3.5) 
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where the quantity F is called the classical Fisher information (CFI). 

The statistical distance depends on the choice of POVM. This dependence is 
to be removed by optimization over all possible quantum measurements. There¬ 
fore, the problem of finding quantum analog of classical statistical distance is 
equivalent to the problem of maximizing the Fisher information over all POVMs, 
i.e., symbolically 

ds%(e(0)) = d9 1 2 maxF(g(6),flA . (3.6) 

fn«} 

The subscript Q reminds one that this is a metric on a space of density matrices 
of quantum states. 

Derived in [39] the upper bound of Fisher information, called the quantum 
Fisher information (QFI) is given by 


-pQ(i?out($)) — 2 


|(i out |5^out(0)|i out )| 5 


^out ^out 
ri ' y j 


= E 


(%>rT 


/„° u t _ „out\2 
2 y- Wi Pj ) i 


^out ^out 


%\ f 


> F(g out (e),fi^) 


(3.7) 


where |* out ) are the eigenstates of g 0 ut (0) with corresponding eigenvalues p° ut 
and the hermitian operator h generates the infinitesimal unitary basis transfor¬ 
mation 


I i out (e + dB)) = e ldeh \i out (e)) = |* out (0)) + ideh\i out (d )). (3.8) 


When QFI is supplemented with an appropriate interpretation and physical 
context it becomes a powerful tool for investigating the structure of particle 
entanglement, as it will be shown in the upcoming chapters. 

The general formula for Fq is quite intimidating and can be difficult to 
work with. Fortunately, this is not always the case. For example, when the 
interaction of the prob and the measured system results in unitary evolution, 
f?out(0) = e~ l9h gi n e l6h , the eigenvalues of input density matrix p‘ n remain un¬ 
altered. Corresponding output eigenstates are related to input eigenstates via 
|i out ((9)) = e~ l6h \i m ). In this simple case the QFI is given by 

Fq{ e~ i0h g, n e ieh ) = 2 ]T \ (H% in ) |' . (3.9) 

P) + Pf 

1,0 L j 

This expression simplifies even further when the input state is pureQ. g- ln = 

l*X*l 


F Q (e~ ieh |W»=4 (T|/i 2 |T) - ((T|/i|^)) = 4(A 2 h)y . 


(3.10) 


Hence for pure states, QFI is simply proportional to the variance of the trans¬ 
formation generator h. 

1 Note that the sum over eigenvalues in Eq. (13.71) also includes cases when one of pi-s is 


zero. 
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I conclude the discussion on distinguishability metric by noting that upper 
bound (1X71) is achievable, i.e. for given g(9) there always exists an optimal 
POVM {fl‘ opt) }, so that F(q(0), f[^ opt ^) = Fq(q(9)) [59]. Therefore, the distin¬ 
guishability metric (13.61) on density operators becomes 

ds 2 Q = F Q (g(9))d0 2 . (3.11) 

The density matrix metric (13.111) also appears in another context. For example, 
a distance between density operators was defined in SB HI S3]: 


d(Qi, Q 2) 2 — 2 


1 — TV 


\J \fdi Q2 \/l >1 


(3.12) 


In quantum information theory this quantity is interpreted as a fidelity - the 
measure of the “closeness” of two quantum states. It can be shown that for 
neighboring density matrices it reduces to 

d(g, g + dg) 2 = ^dsq + o(dg 3 ). (3.13) 


3.3 Precision of estimation 


As I noted above, the problem of precise estimation of 9 is equivalent to the 
problem of distinguishing density matrices along the trajectory g(9). Consider 
now the following procedure: a series of quantum measurements IIj (not nec¬ 
essary the optimal ones) are repeated v times, yielding the results £ 1 ,... 
These results are then used for estimation the value of parameter 9 via a func¬ 
tion called estimator. 9 est = 9 e s t(£i,.. ■, £j,). The variance of any unbiased 
estimator, i.e., functions which satisfy condition 9 est = 9 for all values of 9 (the 
average (...) is taken with probability distributions p(£i|0) ) is bounded from 
below by the so-called Cramer-Rao lower bound [44]: 


A 2 £L 


^ -- 


* F(e 0 ut(e),Ut) 


(3.14) 


According to Fisher theorem S3, for given probability distribution p(£|0) and 
sufficiently large number of repetitions v this bound is achievable by unbiased 
maximum-likelihood estimator |46l l47l . 

The bound (13.141) is known from the classical estimation theory. It can now 
be combined with the quantum result of optimization over POVMs (13.71) . This 
brings us to the final conclusion that the error of estimation of parameter 9 for 
an arbitrary choice of estimator and quantum measurements performed on the 
output state of the probe is bounded by QFI 


A 2 9 > 


1 1 
" F Q (g out (0)) 


1 1 

v dsq/d9 2 


(3.15) 


Aside of its utility for atomic interferometry and quantum metrology as a 
whole, the bound (13.151) also leads to formulation of a generalization of uncer¬ 
tainty principle. Note that QFI can be bounded by a following expression (I 
















26 


CHAPTER 3. EFFICIENCY OF INTERFEROMETER 


drop subscript “out” for convenience) 



(3.16) 


(3.17) 


Here h was introduced in (13.811 and ( A 2 h)$ = Tr (h 2 g) — (Tr (hg)) 2 
on state g. When this bound is combined with (13.1511 we obtain 


is a variance 



(3.18) 


The above inequality transition between a “classical” uncertainty principle, 
when h = 0, which limits distinguishability of probability distributions, and 
a “quantum” uncertainty principle, when dgpi = 0 for all i, which involves the 
generator h “conjugated” to phase 9. 

In summary, I showed that the precision of parameter estimation is tied to 
QFI which quantifies the susceptibility of the probe state to change induced by 
the interaction with measured system. This result is a starting point for the 
upcoming chapters where I describe how the efficiency of metrological task can 
be used for investigation of properties of quantum states. 








Chapter 4 


Dynamical entanglement 


In this chapter I discuss the role of quantum Fisher information as a criterion 
for useful particle entanglement and introduce the concept of dynamical entan¬ 
glement. 

4.1 Quantum Fisher information as an efficiency 
criterion for entanglement detection 

In this section I aim to show that quantum Fisher information (QFI) associated 
with coherent transformation , i.e. transformation which acts on each party in 
the same way, can serve as an entanglement criterion |48j . Moreover, I will 
demonstrate that the value of QFI strongly relies upon the correlations between 
the parties forming the system, thus it can also be used as an indicator of the 
degree of non-classical correlations. 

Consider a iV-parties system that has been initialized in a separable state 
(as in Eq. 


Pin — ((separable — ^ ^ Pk ® ‘ ‘ ® 


(4.1) 


k 


where 0 ^ Pk ^ 1 and Pk = 1- I assume that parties are copies of the same 
subsystem described by a Hilbert space of finite dimension. By doing so I can 
consider a scenario where parties are added or removed from the system. In 
such a case N can be treated as a resource utilized for performing a given task. 

The system undergoes an interferometric transformation generated by co¬ 
herent Hamiltonian H = J2n=i where hf n ' > = <gj /j ® are 

operators acting only in the subspace of party n. The transformation imprints 
information about the phase 9 onto the initial state 


Pin -t f? ou t(0) = e ieH g\ n e x6H 


(4.2) 


Quantum Fisher information associated with this process is given by (see Eq. EB) 



(4.3) 
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where |* out ) and p° ut are eigenstates and corresponding eigenvalues of g out (note 
that the sum also includes cases when one of the eigenvalues is zero). The second 
argument of Fq reminds us that the transformation of input state is generated 
by the Hamiltonian H. I proceed by finding the upper bound on the value 
°f Fq (((separable, H) in order to establish a relation of type (11.111) which is a 
necessary condition for the QFI to be considered as a criterion for entanglement 
detection. 

In section showed that the quantum Fisher information is obtained by 
optimizing classical Fisher information (CFI) over all possible measurements. 
Consequently QFI inherits certain properties of CFI, including convexity as a 
function of density matrix. Therefore, when the input state is separable, i.e. 
it is a convex combination of product density matrices, as in Eq. (ED, i get a 
following inequality 


FQ (Separable, H) < E P k F Q (gj^ ® ■ ■ • ® of\H) 


= 2 E P *E 


|(ig ut |aa(e-** g gj^ ® ... ® g]^V 0H )|jg ut )| S 

P k i+P k j 


(4.4) 


Note that the unitary interferometric transformation does not change proba¬ 
bilities Pfc and eigenvalues pf of product density matrices . The derivative of 
the output state with respect to parameter 6 is given by dg(e~ l0H ge l9H ) = 
—i[H, e~ 10 ^ ge 10 ^ 1 ]. Substituting this result into Eq. (14.41) I get 


k i,j 

= 2 E P *E 

fc 


|(^ ut |[I?, e-^4 1} ® • • • ® Qk N)eie6 ]\jk Ut )\ 2 

Pi + Pj 

l , V m ut \H\f k ut )\ 2 < 

Pi + p* 


< 2 E P *E (p"+p*) m ut \m ut )\ 2 

k i=£j 


4 E Pfc Ep^ A2 ^)*r 


(4.5) 


where (A 2 JJ)jout = (j£ ut | H 2 |i£ ut ) — ((i£ ut | H |*£ ut )) 2 is a variance of H on state 
|i)! ut ). However, unitary transformation generated by coherent Hamiltonian do 
not entangle parties and state g^ ® ® gE remains separable. It follows 

that its eigenstates are product states of form |i£ ut ) = ® ® | <t>^) and 

the variance of Hamiltonian H breaks up into sum of variances of single-party 
Hamiltonians (A 2 H) i out = Y^n= ^A 2 ^)^). The variance itself is bounded by 

the difference of extreme values of operator spectrum, namely (A 2 h) ^ (A max — 
A m in) 2 /4. The bound on variance of single-party Hamiltonian together with 
the conditions YliPi = 1 an<4 Hk -Pfc = 1 yields the inequality which ends the 
inquiry 

Fq (^separable; H) ^ A^(A max A m i n ) 2 . (4.6) 

I proceed to show that non-separable states can give QFI greater than the bound 
above. 

Suppose that the input state is now arbitrary. I start with Eq. m and 
note that this expression can be bounded by the variance calculated with the 
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whole density matrix 


F (h■ - 2 K z ° Ut |^g° u t|j OUt )| 2 _ o ( p ? Ut P°J Ut ) I /'-out| tji 7 '°ut\ 1 2 
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< 2E (p? ut +pD l(* out |^lj out }| 2 = 4^ P r t (A 2 H) i out = 




4Tr (H 2 P° Ut N° Ut ) (*° ut |) - 4 ^ P? Ut Tr (tf|*° ut ) (*° ut |)" 


< 4Tr {H 2 g out ) ~ 4Tr (Hg ont ) 2 = 4(A 2 H% out 


(4.7) 


In turn, the variance is bounded by the difference of maximal and minimal eigen¬ 
values of H which are simply N A max and NX m i n . Thus I obtain the ultimate 
bound on the QFI which cannot be surpassed by any input state 

F Q {Qi n, H) < N 2 ( A max - A m i n ) 2 (4.8) 


Note the difference between bounds 033 and (0H): separable states cannot 
exceed QFI that is proportional to number of parties N while the ultimate 
bound scales with N 2 - the square of number of parties. This shows that indeed 
there might exist some states which give QFI greater then separable bound. 

As a last step, I identify the family of states which saturate the ultimate 
bound to be of “Schrodinger cat” type 


|Cat v ) 


|Amax)® JV +e^|A min )® Ar 

V2 


(4.9) 


where ip is an arbitrary phase and | A max / m i n ) are the eigenstates of single-party 
Hamiltonians h . Indeed, the QFI for such a pure state yields 

FgdCat^H) =4(A 2 H) Ga ^ = 4(Cat v |77 2 |Cat v ) - 4(Cat v |77|Cat v ) 2 = 

= 2 (N A max + N A min ) — (NA max + AA m in) = 

= A 2 ALx + A 2 Amin - 2 A^ 2 Amax A m in = N 2 ( A m ax - A m in) 2 (4.10) 


This concludes the proof but also is a pleasing result in itself. Schrodinger cat 
states, which in case of systems composed of qubits are also called NOON states 
or GHZ states for N = 3, are considered as maximally entangled states (see 
Ref. [49] for extensive discussion on the problem of multipartite entanglement 
measures). Therefore, since maximally entangled state such as cat states can 
yield a tremendous gain of order N over separable states, QFI is also a good 
indicator for a degree of non-classical correlations between parties. 

The relation established by the Cramer-Rao lower bound (see Eq. (13.141) ') 
asserts that QFI is a well defined efficiency criterion associated with task of 
parameter estimation in interferometric transformation generated by H. In this 
context the efficiency introduced in Sec. 11.51 is the inverse of the error of the 
estimation divided by the number of experiment repetitions v 


^estimation of 0)^ = y/ Fq (g. m , H) . 


(4.11) 
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For a separable state the efficiency is no greater then y/N |A max — A m i n |, or the 
error scales as A 0 ~ 1/y/N - the precision reaches the shot noise limit. Cer¬ 
tain entangled states allow for beating this limit and for maximally entangled 
Schroedinger cat states the improvement is of order Ad ~ 1 /N, which is par¬ 
ticularly lucrative for large number of parties. The ultimate bound (14.81) which 
defines the maximal achievable precision is called the Heisenberg limit (HL). 

Quantum Fisher information backed up with the context of atomic inter¬ 
ferometry will now be the centerpiece of the upcoming discussions. I will use 
this powerful concept to analyze and classify non-classical correlations present 
in ultra-cold atom systems. 


4.2 Introducing dynamical entanglement 


The troublesome aspect of using the QFI as a measure of entanglement is the 
requirement for associating it with some transformation which imprints param¬ 
eter to be estimated. Of course, the context of atomic interferometry provides 
us with the choice of transformation. Nevertheless, still it would be desirable to 
make an attempt to find some generic type (or types) of transformation which 
would justify the use of QFI on its own. 

In chapter [2] I argued that treating atoms as qubits provides a satisfactory 
description of ultra-cold atom systems used for ongoing research on quantum 
entanglement. By restricting the discussion to qubit system I will be able to 
detach QFI from the context of interferometric experiment by removing the 
ambiguity of choice for coherent transformation. Then, the QFI becomes a 
function of the state only and it can be treated as a quantity characterizing its 
properties, including non-classical correlations. However, the relation between 
QFI and atomic interferometry is too valuable and it would be unwise to discard 
it all together. As I will argue, the direct correspondence between this generic 
QFI and arbitrary interferometric transformation can be easily restored. 

From now on I shall assume that parties constituting the system of interest 
are particles living in a two dimensional Hilbert spaces. This restriction sim¬ 
plify the problem significantly. In case of qubit systems single-particle transfor¬ 
mations, which build up coherent transformations, are mapped onto rotations 
of spin 1/2 system generated by triple of spin operators J- n ^ = ct^/ 2, with 
i = x,y,z and a- TI> are Pauli matrices of n-th qubit. Any given unitary trans¬ 
formation of a qubit can be characterized by an angle 9 and the unit vector 
n parallel to the axis of rotation in a following way: e~ l6n J . However, the 
reference frame can be changed at will and it can always be chosen in such a 
way that the new z-axis coincides with the axis of rotation n. Therefore, each 
transformation can now be considered as a family of rotations about the z axis 
characterized by a single parameter - the angle 6. Now imagine that such a 
generic transformation describes an interferometer while the angle 6 is an un¬ 
known parameter to be estimated. The quantum Fisher information associated 
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with the precision of this fictional estimation task is given by 


Fq(Qi n, Jz) = 2 y 


|(i OUt |[^,eo U t]|j° Ut }| S 


= 2E 


h3 




iPi -PjY 

Pi + Pj 


KW)I 2 = 


= 4^2pi{A 2 J z }i - 8 ^ 




PiPj 
Pi + Pj 


\m\j)\ 2 , 


(4.12) 


where g out = e~ l9Jz g- 1 ^e ieJz and p° ut are its eigenvalues with corresponding 
eigenstates |i out ). In order to obtain this equation I used the fact that unitary 
transformation does not change eigenvalues of density matrix ( p° ut = pi) and 
(* out |J s |j out ) = (i\e l9 ^ z J z e~ l9 ^ z \j) = (i\J z \j), where |z) is an eigenstate of />j n 
corresponding to pi. As a result, Fq does not depend on the value of 6 and 
the explicit form of the transformation e~ l9Jz becomes superfluous. Therefore, 
I introduce new function of state inspired by QFI, which stands on its own, not 
relaying on a context of interferometric experiment 

He) = 4^P.(A 2 i z ) i - 8 E ^Srl<*l^l-?')l 2 • ( 4 - 13 ) 

i& Pi+Pj 

I shall refer to F as dynamical susceptibility. The name is inspired by the 
interpretation of QFI as a susceptibility of the state to change due to driving 
by the system Hamiltonian, as it was discussed in Ch. [3] 

Since T equals QFI associated with interferometer generated by e~ l9Jz all 
the properties derived in previous section carry over. Hence the dynamical 
susceptibility can also serve as a criterion for entanglement detection. Since 
spin operator a z /2 has only two eigenvalues +1/2 and —1/2 a general SNL 
bound (14.611 reduces to 

-F(^separable) A^ . (4.14) 

The ultimate HL bound (14.811 in this case is 

F{g) < N 2 . (4.15) 


Similarly to criterion based on QFI, dynamical susceptibility will serve as the 
mean to analyze and quantify non-classical correlations in ultra-cold atoms sys¬ 
tems for which qubit approximation is valid. I will say that a state g which 
gives F(g) > N is dynamically entangled to the degree F(g). 


4.3 Relation between useful and dynamical en¬ 
tanglement 

The aim of this section is to establish a procedure to relate the dynamical 
susceptibility with an interferometric experiment. To achieve this goal I need 
to be able to express the QFI associated with a task of estimating a value of 
phase imprinted by an arbitrary interferometric transformation in terms of F. 

A general (unitary) interferometric transformation can be parametrized by 
a sequence of rotations with one of the angles being unknown 0 

—. e ~ioi 1-j e -ja fc • jg-ie n- jg-i/31 • j e -iPm ■ 3 _ 

—ioL • J e ~i0 n • J e ~iP • J = fjn ■ 3 fj^ 


= e 


(4.16) 
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Here U a /p is a composition of rotations (therefore being rotation itself) preced¬ 
ing/following the phase imprint. The task is to estimate the value of 9. As 
usual, the precision of estimation is bounded by QFI, which in this case reads 


F Q (g in ,U e ) = 2j2 

i,3 


|(» out |^gont(g)|j out )| 2 _ |(i out |[n' • J, g ou t]|j' out )| 2 

Pi+Pj Pi+Pj 

{Pi -Pj) 2 l/joutm/ . j|jOut)|2 (4.17) 

Pi+Pj 1 1 


Here n' = R Q 1 n is a unit vector n rotated counter-clockwise by angle |a| about 
the axis a by means of rotation matrix R Q . This relation follows frond 

deBont = {deUe)QiJjg + Ugg in (d e Ul) = 

= (-in - (U a 3Ul))g out + g out (i n • (U a 3Ul)) = 

= -i[ (R" 1 !!) • J, g out ] = -i[ n' • J, g out ] . (4.18) 


Further manipulations confirm that Fq is independent of 6 


f q = 2 J2 


c Pi-Pjy 

Pi+Pj 


|(z out |n'.J|j out )| 2 =2£ 


( Pi-Pj ) 2 

'ti Pl+P > 
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{Pl Pj)2 Ki|(R^)-JU)| 2 = 2£ (Pi Pj)2 |<»|n w -J|j)| 2 (4.19) 

Pi+Pj £ Pi+Pj 


Where R/j is the clockwise rotation about axis (3 by an angle |/3| .The final step 
before I can relate Fq to dynamical susceptibility is to introduce a rotation 
R such that R(n” ■ 3)R^ = J z . This transformation always exist and can be 
applied in (14.191) 

F Q (Qin,Ue) =2 £ ~ Pj)2 \(i\&R(n"- 3)&R\j)\ 2 = 

Pi+P > 

= 2j2 iPl ~/ :)2 = 

& Pi+Pj 

= 4£p,(A 2 J z ) m - s£ ^^\((i\tf)j z (R\j))\2 = F(Rg in R t). (4.20) 

, , IP ~ P:i 


Hence, any usefully entangled state can be transformed by means of rotation 
(which does not introduce any entanglement) into state that is dynamically 
entangled to equal degree. The particular transformation is uniquely determined 
by the interferometric sequence in question. 


1 One can consider even more general interferometric transformation, where the phase 9 is 
imprinted simultaneous with another rotation, i.e. 

'9 


Ug = UaTex p 


—i / ds 
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l + 7(s) J 


Up = UaU-yT exp (—i9£ ds %(rs){n ■ J)C/+ (rs)^ Up 


where 7~exp is a time-ordered exponential and = 'T exp ^—2 fg ds'j(s) • Jj . In such a case 
we have dgg ou t = — i[(Ra 1 R£ 1 n) • J, g Q ut], with n- J = f* ds t/ 7 (rs)(n • J)I/^(rs). 
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The reciprocal relation also exist but is not unique. Working out the deriva¬ 
tion backwards I obtain the following relation valid for arbitrary rotation R 

He) = Fq{RqR\ U a U n e- ieS *UtUp). (4.21) 


Rotations U a and Up which precede and follow imprinting of the phase along 
axis set by U n satisfy condition R = U^U^Up. That is, these transformations 
are decomposition of R into sequence of rotations. Such a decomposition can 
be done in infinite number of ways. The only constrains can be imposed by a 
practical concerns laid down by a context of experimental setup. 

Finally, I consider a case of non-unitary interferometric transformation. Such 
transformation not only rotates the eigenstates of density matrix but also ma¬ 
nipulates its eigenvalues which describe classical ignorance of the observer re¬ 
garding the preparation of the state. In principle, the information about the 
parameter to be estimated can be drawn form both of these processes as it 
is demonstrated by the QFI for a general non-unitary transformation which 
depends on the parameter 6 


7-inon—unitary 


! E 


cp° ut - vTf 

it _|_ ^out 
' Pi 


|(* out |h|j out )| 2 + ^ 


(d e p° ut ) 2 

^out 

Pi 


(4.22) 


where the matrix elements of hermitian operator h describing the “unitary” 
part is defined as a generator of the infinitesimal transformation |i out (0 + (i0)) = 
e~ ld0 ^\i ont (6)) (as in Eq (13.81) 1. In general -pQ° n ~ unltary j s no t a p r0 per effi¬ 
ciency criterion since relation such as m does not exist. This is because only 
the “unitary” part depends on non-classical correlations present in the system. 
Therefore, QFI associated with transformations that are not unitary can be 
biased by the additional “classical” part. The dynamical susceptibility is not 
burdened by this flaw, as by its definition, it “picks” only the unitary part of 
the QFI. By choosing R such that the generator h is transformed into J z we 
obtain the following relation 

F non-unitary = t) + £ • ( 4 - 23 ) 

i Pi 


We see that the dynamical susceptibility indeed describes the “unitary” part 
of Fg° n_unltaiy . Note that T no longer depends on the input state g- ln only. 
Instead it is a function of the output state g out which means that it depends on 
the type of transformation, the duration of the experiment as well as the initial 
state. This is not a surprise because in general non-unitary interferometric 
transformations include the process of decoherence which destroys correlations 
within the system. Therefore, the degree of entanglement is no longer conserved 
and this is reflected in the dependence of dynamical susceptibility on the details 
of the process. 

Relation between dynamical susceptibility and the QFI associated with a 
given transformation discussed above can be understood as follows. The inter¬ 
ferometric experiment consists of three stages: state preparation, phase imprint 
and measurement. So far we focused on the phase imprint stage represented 
by the interferometric transformation Ug = U a U n e~ ieJz U} l Up■ We took it for 
granted that the preparation stage already took place and the result was the 
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input state pi n . Also we were never concerned with measurements which follow 
the imprint because QFI is optimized over all possible realization of this stage. 
However, the boundaries between the stages are not clear-cut. Formally one 
can regard the transformations U^U /3 and U a U n which precede and follow the 
phase imprint as a part of state preparation and a measurement stage, then 
e -iQJz becomes the new interferometer. In fact, the relation (14.2011 is an exam¬ 
ple of this formal division. Indeed, dynamical susceptibility is equivalent to the 
QFI associated with a generic interferometer e~ l9Jz and the transformation R 
represents the net result of transformations preceding and following the phase 
imprint. 


4.4 Dynamical entanglement and spin squeezing 

The spin squeezing is another efficiency criterion for entanglement which is re¬ 
lated to precision of two-mode interferometer. Spin squeezing parameter is 
defined as the ratio of the error of the phase estimator derived from the popu¬ 
lation imbalance measurement in the Mach-Zehnder interferometer to the SNL 
precision. Before I compare criterion based on QFI and the spin squeezing I will 
review the derivation of £s- 

In the language of pseudo-spin, the Mach-Zehnder interferometer is repre¬ 
sented by the following sequence of rotations 

i/jMZ) = Jx (4.24) 

and as usual 9 is an unknown phase to be estimated. The rotations about x-axis 
which precede and follow the phase imprint represent beam-splitter operations 
which enable “mixing” of atoms occupying the modes of interferometer. The 
value of 9 can now be inferred from the oscillatory behavior of population im¬ 
balance between the two modes of the output state g 0ut = Ug Z ’g- m Ug Z ”. 
According to what I showed in Sec. 12.31 the 2 component of the spin operator 
is an observable associated with the difference of mode populations. We can 
easily verify that the expectation value of J z indeed depends on the parameter 
9 


{Jz) 


out 


Tr[j z (U^g in U^)\ = Tr [(U^ J z U^)g in ] 
COS0(J,)in + sin0( J x )in = f (9) 


(4.25) 


The estimator for parameter 9 can be chosen so that 9 es t = / _1 ((</ 2 ) out ). The 
uncertainty of this estimator is determined by the variance of the population 
imbalance and it can be calculated using the error propagation formula: 
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(4.26) 
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where v is a number of experiment repetitions. The spin squeezing parameter 
is defined as 


£s(e in) = v 


( lim e ^o A 2 0 

V v n 



, r (A 2 i,)in 

I < Jx } in | 2 


(4.27) 


Spin squeezing is also a criterion for particle entanglement mum os]: 


£,S (^separable) 
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ZUe) 




(4.28) 


The main advantage of the spin squeezing is how relatively easy it is to as¬ 
sess in most experimental settings. The variance of population imbalance 
((A 2 J Z )) can be deduced from standard particle number measurements which 
are always setup so that they can resolve between different modes of the inter¬ 
ferometer. The denominator, |( J x )\ measures the coherence between two modes. 
It can be obtained by examining the visibility of interference fringes observed 
after the “mixing” of matter waves from each mode. For example, this can be 
achieved by releasing the BEC from the trapping potential and allowing the 
atomic clouds to freely expand and overlap thus forming the interference pat¬ 
tern. High resolution imagining techniques allow for very precise measurements 
of the structure of this pattern [33] ■ In contrast, extracting QFI is also possible 
but it requires highly sophisticated methods |50) . 

The “standard” choice of axes defining £$ came out naturally as a direct 
result of particular estimation strategy adopted for the Mach-Zehnder interfer¬ 
ometer. The first generalization of the spin squeezing parameter is to replace 
the denominator with the square of the length of projection of average spin 
vector onto the plane perpendicular to the 2 axis: 


Zs = 



(A 2 Jz)in 

| ( Jx ) in | 2 4“ | ( Jy ) in | 2 


(4.29) 


A more general definition which explicitly points out the flexibility in rearrang¬ 
ing the directions of squeezing is given by 


Z'sie, n) 



(A 2 (n ■ J))g 

ll(J)±ll 2 


(4.30) 


Here ||(J)x|| denotes the length of vector (J)j_ which is a projection of the av¬ 
erage spin vector (J)g onto plane perpendicular to n. Similarly to relations be¬ 
tween dynamically and usefully entangled states, a state which is spin squeezed 
in respect to one direction ni is related to state squeezed in direction ri 2 through 
unitary rotation 

(A n i) = £ 5 (RgR ?, n 2 ), (4.31) 

where .ff(ni -3)R^ = n 2 J. In particular £s{q) = £g(f?, e z ) = £s(g) if the axes 
of reference frame are chosen so that e^, || (J)_l- These relations are analogical 
to relations I have established between QFI and dynamical susceptibility which 
is essentially the QFI for “standard” interferometer. However, the components 
of pseudo-spin do not refer to direction in real space and J z is in fact the popu¬ 
lation imbalance operator while J x is related to visibility of interference fringes. 
Passing to the rotated frame of reference might sabotage the main advantage 
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of spin squeezing parameter: the ability to measure it in experiment. There¬ 
fore, in almost all circumstances the states which are spin squeezed according 
to “standard” parameter (be- such g that (g) < 1) are most desirable for 
practical uses. 

From theoretical point of view, dynamical susceptibility (and QFI) is more 
attractive then spin squeezing. Since is a variance of particular estimator 
it is bounded from below by the inverse of QFI associated with Mach-Zehnder 
interferometer 


A 2 0 est 


1 ^ (gin) 1 1 

* N " u F Q (g in ,U^ MZ) ) 


F Q (g in ,Ur ] ) 
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1 

if® 


(4.32) 


Therefore, QFI, and by extension dynamical susceptibility, detects more types 
of entangled states then spin squeezing |50l . The on-site atom-atom interaction 
(as it is the case for atomic interferometers realized in double well setup) is the 
natural source of spin squeezed states (i.e. states for which < 1 is satisfied). 
Also this type of states have the greatest overlap with usefully and dynamically 
entangled state detected by QFI or dynamical susceptibility. Nevertheless the 
overlap is not perfect. In [48] authors carry out a detailed analysis of this family 
of states and showed that indeed QFI criterion is able to detect more entangled 
states then spin squeezing. 

The advantage of methods for entanglement classification and analysis based 
on QFI is most explicit when one considers states created in atomic collisions 
(as in the case of twin-beam type of atomic interferometers). The characteristic 
feature of this type of states is the vanishing of the average spin vector (J) 
[35]. This leads to zero visibility of interference fringes which renders the spin 
squeezing parameter undetermined because the denominator in Eq. (14.2711 equals 
zero. Similar problem is encountered in case of Schrodinger cat states. This is 
indeed a serious drawback since these states are considered to be very strongly 
entangled. It seems reasonable to stipulate that the analysis of this type of non- 
classical correlations might be the key to understand the nature of entanglement. 









Chapter 5 

Physical interpretation of 
dynamical entanglement 


In this chapter, I employ the dynamical susceptibility, its interpretation as a 
measure of the susceptibility of state to change, and the equivalence between 
bosonic qubit- and spin-systems to discuss the possible physical interpretation 
of particle entanglement encountered in cold-atom systems. 


5.1 The role of particle indistinguishability 

The key aspect of ultra-cold atom system is the indistinguishability of bosons 
which constitute it. It is of great importance to analyze this contribution to 
overall non-classical correlations present in these systems. To this end I will 
exploit the correspondence between system composed of N qubits and a system 
composed of a single pseudo-particle with spin degree of freedom. Initially 
I will not assume that qubits are identical bosons and I will investigate the 
changes in the degree of dynamical entanglement when the indistinguishability 
is imposed. By doing so I will establish important relation between the spin of 
pseudo-particle, entanglement and the indistinguishability of particles. 

In Sec. EH showed that no state can surpass the Heisenberg limit (HL) of 
dynamical entanglement. I have also found the class of Schrodinger cat states 
which reach this limit. Now I shall examine the role of particle indistinguishabil¬ 
ity in attaining the ultimate HL. I start with the following bound on dynamical 
susceptibility 


m = 4^><A 2 J^ - 8^ 

i 7 £P' + p> 

< • (5- 1 ) 

i i 

With this bound I can utilize the equivalence between the qubit- and the spin- 
system. On one hand, eigenstates | i) can be written as a superposition of 
products of single-qubit basis states, i.e. | i) = Cal...a N \ a i) ® ® |<tjv), where 

a n =t)4~ On the other hand, each qubit is equivalent to spin 1/2 particle and 
the basis states correspond to spin eigenstates: |t / 4-) ^ |s = \,m = ±i). 
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Since states |i) are superpositions of outer products of N spin eigenstates they 
can be expanded in the basis of irreducible representation of the rotation group, 
i.e. the basis of total angular momentum, 

l*> = li m Ar) = 

— J — 

= EE E = (5.2) 

J =0 v M——J J=0 

States | are normalized projections of |*) on the subspace with the total 
angular momentum J , while v is a degeneracy index labeling representations 
with the same J. Note that in general the sum over total angular momenta J 
ranges from 0 to N/2. 

Subspaces with definite J are invariant under rotations, thus the expectation 
value of J 2 on the state |i) is equal to the sum of expectation values on each 
projection i.e. (i\J 2 \i) = J2j\Cj^\ 2 (ipj^\J 2 \iJjj^). In the final step, I 

note that (i/)j^\j 2 \'i/)j^) ^ J 2 , so the bound for dynamical susceptibility from 
Eq. (15.11) is 


N_ JV 

,F(£K4;£J 2 5> \C^\ 2 =4j2 j2 Tr(iljg) . (5.3) 

J=0 i J= 0 

Here flj = \ {J,M-,u\ are projectors onto subspace of total 

angular momentum J. This bound is more strict then the ultimate HL it 
takes into account the symmetry properties of the state through probabilities 
of finding the system in subspaces of total angular momentum 

N_ 

2 

Pj = Tr (Ujg) => Pj ^ 0 and ^ Pj = 1. (5.4) 

j=o 

Bound (15.31) is an average square of total angular momentum J 2 weighted by 
probabilities Pj. If the state of the system is spread among wide range of to¬ 
tal angular momentum subspaces it is less entangled then the state which is 
distributed only among small number of subspaces with the highest J-s. In 
extreme case, the entanglement of a state confined to J = N/2 subspace is 
bounded by the HL of N 2 . Moreover, the projectors II j are rotationally invari¬ 
ant, i.e. for any rotation R we have [IIj,i?] = 0 and hence Pj = Tr(IIjp) = 
Tr(IIjI?^I?p) = Tr(IIj(.RpiH)). Therefore, the efficiency of the parameter esti¬ 
mation in interferometric experiment is also bounded by the decomposition of 
the state into total angular momentum subspaces since according to Eq. (14.201) 

^ ^ N_ 

we have Fq(q, U e ) = P(RgR^) < 4 J2j=o p J j2 - 

The sum over J in Eq. results from the rules of addition of angular 

momenta of the pseudo-spins of qubits constituting the system. Now I will 
review how the total angular momentum eigenstates \J,M\v) are constructed 
out of product states of qubits’ pseudo-spin eigenstates. By examining this 
process we shall gain a necessary insight to answer the question what role is 
played by the indistinguishability of qubits. 




5.1. THE ROLE OF PARTICLE INDISTING UISHABILITY 


39 


I start with the even set of N qubits (the case of odd N is solved analogically). 
I take a symmetric superposition of first J—M qubits in state 11, — |) and J+M 
qubits in state ||,|) then multiply it by -y — J anitsymmetrized pairs of qubits 
that remain. Such a state is a total angular momentum eigenstate: 

\J,M-u= 1) =5 . (5.5) 


Here S is a symmetrization operator and 


l*-> 


\b\) 


\ii _i\ _ 11 _i\ii i\ 
/I 2’ 2/ I 2’ 2/I 2 5 2 / 


V2 


(5.6) 


is the antisymmetric singlet state. The value assigned to the degeneracy index 
v indicates that this particular choice of which qubits are to be antisymmetrized 
is one of many possible^ From this construction we see that the total angular 
momentum is set by the number of symmetrized qubits forming the state. In¬ 
deed, when a total angular momentum operator, + = Y^,n=i '^k ' > - ac ^ s on this 
state, the singlet pairs of qubits never contribute: 

J k (S = 

+ S[ E (5-7) 

\n=2J+l J 


and now I can rearrange terms in the sum of the second line of the above 
equation 
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The last equality follows from the fact that singlet state vanishes under action 
of any global angular momentum operator: 


(ji") + jj" +i ))+_} = (+±-±) 

) li± 


lid) 


V 2 


= 0, 


(j£ n) + i£ n+1) 


= 0. 


(5.9) 


Hence, all that remains from Eq. E3 is the action of the operator on the 

x It can be show that the number of distinguished choices equals d J< N_ = (n/^— j) ~ 
( N /2-J-l) atld d f = L 
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symmetrized part. An elementary calculations show that 


n =1 
2 J 

n=l 

= y/(J ± M)(J ± M - 1) 5 [IA, 1 )®^+M±ij 5 (5. 10 ) 

which confirms that Eq. (1531) is correct. 

Recall that the state describing a system of identical bosons has to be sym¬ 
metric in respect to qubit permutations. Suppose that N' < N qubits are 
indistinguishable while the remaining N — N' particles remain distinguishable. 
In that case the minimal value of the total angular momentum that can be at¬ 
tained by the system is restricted to ,/ ^ 4J-, because any possible state has to 
be symmetrized in respect to at least N' qubits. Therefore, when all N qubits 
are identical, the bound (15.51) grows to become the Heisenberg limit J-{g) < N 2 
since only maximal J = N/2 is allowed, which implies that Pj—n/2 = 1 an d 
all other Pj are zero. Hence, the indistinguishability of qubits “automatically” 
confines all possible states into subspace of maximal total angular momentum 
which enables potentially highest degree of entanglement. If particles where not 
identical, it would require a tremendous amount of effort and ingenuity to be 
able to perform an experiment with a large ensemble of qubits prepared in a 
symmetrized state. For example, Schrodinger cat states (see also Eq. (14.91) 1 
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(5.11) 


yield HL degree of entanglement and are explicitly symmetric in respect to qubit 
permutations even if the particles are distinguishable. However, it is well known 
that in practice it is very difficult to prepare such a state for a large number of 
particles. 

The most prominent example of “experimental friendly” state is a so called 
twin Fock state m- In two mode approximation such a state is often described 
as a separable in modes ket (here L/R indicate the two modes of the 

system). In the language of pseudo-spin it is given by a symmetric total angular 
momentum eigenstate 


. N 


\N\ \N\ _ cHl _1\®'T|1 1\< 

I 2 ' L\ 2 /R U 2 ’ 2/ 2 I 2 ’ 2 / 


^]=lf,0>. 


(5.12) 


Such a state is relatively easy to obtain in various experimental setups. For 
example, the twin Fock state is obtained when the potential trapping Bose- 
Einstein condensate is adiabatically brought into a double well trap |34| . Due 
to repulsive interaction between bosons, the condensate is split evenly between 
the wells. 

Although separable in modes, the twin Fock state is usefully entangled to 
a very high degree since QFI associated with Mach-Zehnder interferometer 
t /g MZ ) = for this state is of order of half HL. This interferom¬ 
eter utilizes the beam-splitter operation to prepare the state, therefore 
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the state corresponding to Twin Fock which posses equal degree of dynamical 
entanglement is |TF) = 


J'(|TF))= J F Q (|f,0),17 e (MZ) ) =jv(y + l) . (5.13) 


The entanglement between particles forming this state can be ascribed solely to 
indistinguishability of qubits. Indeed, consider an experiment where instead of 
adiabatic splitting the double well setup is created by bringing together inde¬ 
pendently prepared condensates. If the number of particles in both wells is ex¬ 
actly the same, say N/ 2 , and the two-mode approximation holds, then the state 
of the system is again |TF). This result might be counter-intuitive, since one 
might expect a separable state because the contents of each well never interacted 
with each other. However, the particles forming both condensates are identical 
bosons and the state has to be symmetric in respect to qubit permutations, 
therefore it cannot be separable. To see this lets denote by |f / 4-) the state of a 
particle being localized in left/right well. If particles were not identical the state 
of the whole system would be a product |t) 2 | 4 -) 2 , but for identical bosons 

the state is symmetrizecH | = <S[|^, — 5 )®^ |^, 5 )®^] = 1 ^, 0 ) = 

e -jf J* |TF). I again underline that the beam-splitter e l ^ Jx is a coherent trans¬ 
formation and it does not introduce any entanglement between particles. The 
only reason why dynamical susceptibility surpasses shot noise limit is due to 
symmetrization enforced by indistinguishability of bosons. 

The physical interpretation of this seemingly unnerving result is as follows. 
Although each condensate is prepared far apart from each other, say one was 
created on the Moon while the other remained on Earth, it is incorrect to assume 
that identical bosons forming them where not correlated. Since these particles 
are indistinguishable it is impossible to tell which atoms where on the Moon 
and which on the Earth, therefore the correlation always exists. It follows 
that the symmetrization is enforced even if the two condensates where never 
brought together. Now the question is whether the system of Moon and Earth 
bound condensates is entangled. Formally one can define rotation operator 
gif Ja, representing the beam-splitter and calculate dynamical susceptibility for 
this Moon-Earth double well system. However, it is not correct to claim that the 
state is dynamically entangled. The distinction between artificial entanglement 
such as this, and the proper, physically meaningful entanglement can be made by 

2 In realistic circumstances it is impossible to predict how many particles will condense 
in a given run of the experiment. As a result the state of the whole system is a statistical 
mixture of form 

fco = x) [it) 0JVt i4)®^]5[(ti®^ ai®^] 


Here p(N) is the probability of creating a condensate of N particles. Here I assume that the 
probability distribution p is flat on interval [AT m i n , 7V max ]. The dynamical susceptibility and 
QFI of this state are 


'-Stwoe-** ■ } *) = F 0 (etwc,f/r Z) ) = N (f + l) [l-(^) 2 ] - i(s^5) 2 > N 


Here N = Nmax + Nmin and AN = TVmax — N m i n . Since T (and Fq) surpass SNL the state 
is almost always entangled (with the exception of maximal AN = 2N when T = N). For 
example, when AN = 25%TV and N = 1000 we have T ~ %N. 
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invoking the context of interferometry. Dynamical susceptibility for this system 
is meaningless since it is impossible to create an interferometer operating on 
remote condensates and the relation such as (15.131) between QFI and T cannot 
be established. Therefore, it is necessary to bring the condensates together and 
enable them to “mix” with each other, otherwise the non-classical correlations 
due to indistinguishability are impossible to observed and to utilize. 


5.2 Symmetry and the dynamical entanglement 


In the previous section I have established a relation between the total angu¬ 
lar momentum of a multi-qubit system and the maximal degree of dynamical 
entanglement this system can achieve (see Eq. (15.31) 1. I have found that only 
completely symmetric states, i.e. states with maximal total angular momen¬ 
tum, are able to reach the ultimate Heisenberg limit of entanglement. The key 
question is which properties of the state allow for tapping into this potential 
entanglement that can be detected by a criterion such as the dynamical suscepti¬ 
bility. Once I identify this property, I will be able to make a positive statement: 
we will know what entanglement is , not just what it is not. 

I begin by revisiting the example of Twin Fock state, which written in the 
basis of pseudo-spin eigenstates is given by 


|TF) = e^ J *|f ,0) = C 0 |f ,0) + £ C, 


I N 


2 ) 2m ) 


I N 


2 . ~ 2m ) 


m —1 


V2 


(5.14) 


Here C M = V2D& 0 (-%,- f,f) is an element of Wigner D-matrix[f] We see 
that the state |TF) is a superposition of “Schroedinger kitten” states of form 
(|f, M) + |f, — M))/y/2. Similarly to Schrodinger cat states, the dynamical 
susceptibility of the kitten state is determined by the difference of angular mo¬ 
mentum projections of eigenstates forming the superposition, 


/|f,M) + |f,-M)\ 

l ^2 ) 


4M 2 . 


(5.17) 


The dynamical susceptibility of superposition of kitten states (15.141) is 

2 

J-(|TF)) = 4^ M 2 |C m | 2 - (5.18) 

M—0 

A simple intuitive picture can be draw to interpret this result. The density ma¬ 
trix of a kitten state when written in the basis of angular momentum eigenstates 

3 A rotation operator can be written as 

R = e- iaJx e-^ J ye~ i ' yJz , (5.15) 

where (a;, /3,7) are Euler angles. The Wigner D-matrix is a square matrix of dimension 2 j + 1 
with general element 


(5.16) 
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is of form 


V2 V2 

( '. 


1 

2 


V 


0 0 
0 1 0 
0 0 


0 0 
0 1 0 
0 0 


\ 

0 0 

010 ... 

0 0 


0 0 

0 10... 
0 0 


(5.19) 


The only non-zero elements are found on the diagonal and in the band located 
M units away from the diagonal. The dynamical susceptibility can now be 
broken down as follows: T = 4(2 x M 2 x 4). It is equal four times square 
of distance from diagonal to upper matrix element (M) times the magnitude 
of this element (1/2) plus the same for lower matrix element. This formula is 
analogical to the classical moment of inertia of two points of mass 1 /2 positioned 
at distances M away from the rotation axis. This analogy can even be extended 
to the state |TF), since Eq. (15.1811 resembles a formula for the moment of inertia 
of a collection of pairs of points corresponding to kitten states with masses 
given by probabilities |Cm| 2 /2 and the distances given by the projections M. 
Although the elegant form of Eq. (15.181) results from the a particular form of |TF) 
being a pure superposition of kitten states, the analogy between the moment of 
inertia of density matrix elements and the dynamical susceptibility is generally 
valid. Below I show how it comes by and for that I shall start by introducing a 
language of spherical tensor operators. This language enables a convenient way 
for quantification of the symmetry properties of a state which, as it will turn 
out, is the main ingredient of the dynamical entanglement. 


5.2.1 Basis of spherical tensor operators 

In order to properly describe the rotational symmetry of a state its density 
matrix has to be decomposed into parts which transform independently under 
rotations. For example, recall that the density matrix of a single qubit, which is 
equivalent to spin 1 /2 system, can be written as a combination of spin operators 
and the unit matrix: 

^single qubit — 1 "b 2(J) * «J , (5.20) 

where Ji = \<7i and (Jf) = TV ^ Jiftsingie qubit) • An arbitrary rotation can only 
transform spin operator into combination of spin operators, while the unit ma¬ 
trix is always unaffected. Therefore, (15.20[) is the desired decomposition, where 
the part spanned by triple of J^-s is independent of part spanned by 1. Although 
the decomposition (15.201) might seem to trivially result from the properties of 
Pauli matrices cq, in fact it is a consequence of a more general principle. In 
the language of group theory the spin operators are vectorial objects. To be 
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more precise, they are proportional to spherical tensor operators [511 1521 153] 
with angular momentum j = 1 and magnetic number m taking on values of — 1, 
0 and 1: 


Tj= i, TO =o — Jz (5-21) 

Tj=l t m=±l = (j x ± iJy) = =F;^=«4 (5-22) 

The spherical tensor operators Tj m are an operator analog of the angular mo¬ 
mentum eigenstates \j. m). For example, one way of defining the spherical tensor 
with quantum numbers j and m is to require for it to satisfy a following set of 
equations 

[Jz,Tjm\ = mTjm (5.23) 

[J±,Tj m ] = \J ( j ±m)(j±m+ l)f jm ±i , (5.24) 

which imply 

E [J i ,[Ji,T j m}}=j(j + l)f jm . (5.25) 

i=x,y,z 

These equations are in the direct correspondence with a similar set satisfied by 
the angular momentum eigenstates: 

Jz \j, rn) = m \j, m) , (5.26) 

J± | J, m) = \J ( j ± m)(j ±m + l)\j,m±l) , (5.27) 

J 2 \j,m) =j(i + l)|i,m). (5.28) 

Since the action of an arbitrary rotation on an operator is determined by com¬ 
mutators such as (15.231) and (15.241) it follows that the transformation properties 
of spherical tensors and angular momentum eigenstates are the same. One of 
the consequences is that any rotation transforms a spherical tensor with angular 
momentum j into a combination of spherical tensors with the same j. 

The set of spherical tensor operators forms an orthogonal basis with respect 
to the scalar product (A\B) = Tr(A^B) in the space of density matrices of 
systems with fixed angular momentum. In particular, the density matrix gj 
of system with angular momentum J, operating in 2J + 1 dimensional Hilbert 
space, can be decomposed into linear combination of spherical tensor operators 
with angular momenta ranging from j = 0 to j = 2 J: 

2 J j 

QJ = ^2 'Yl OrnTjm , (5.29) 

j=o m=—j 

Here Tj m = (T^ m |Tj- m ) -1 (Tj- m |gj) are the expectation values of the spherical 
tensor operators. The parts of gj spanned by spherical tensors with fixed j 
have different symmetry and they transform independently. Therefore, decom¬ 
position such as this provides a proper quantification of symmetry properties 
of a state described by gj. Since the density matrix of N bosonic qubits is 
equivalent to a matrix of a system with total angular momentum J = N/ 2, this 
decomposition provides an excellent tool for carrying out my analysis. 
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5.2.2 Quantum moment of inertia 

In order to formalize the analogy with the moment of inertia I shall utilize the 
decomposition (15.2911 in conjunction with a theorem which relates dynamical 
susceptibility with a variance of spin operators. 

For an arbitrary pure-state ensemble {pk,\4>k)} of a density matrix g (i.e. 
0 = X] fc Pfc|0k)(0fc| an d Pk ^ o, J2kPk = 1) the dynamical susceptibility is 
bounded by the average variance [Mi '55] : 

^(p)<4^p fe (A 2 J,) 0fc . (5.30) 

k 

We have already seen an example of this relation when the ensemble in ques¬ 
tion was composed of orthogonal eigenstates of g, for example in Eq. m - 
The theorem states that among all possible ensembles there exists a special 
one which minimizes the average variance which coincides with the dynamical 
susceptibility [55], therefore 

HO) = , m in 4^p fe (A 2 J 2 ) 0fc . (5.31) 


Note that the variance of j z on any pure state can be written in terms of 
expectation value of the commutator with the density matrix of the state 

He) = - min 4^]p fe (A 2 J z ) 0fe = 

= min 2 ^p fe Tr (|0 fc )(^ fc |[j z , [J z ,|^ fe )(0 fe |]]) . (5.32) 

Each pure-state density matrix can now be decomposed according to Eq. (15.2911 , 
J2 3 m=-j T jm ) Tjm- Using the fact that spherical tensors are 
eigenoperators of J z [see Eq. (15.2311 ] we get 


He) 


N j 

T u 2 E ft E E ) ^(\<t ) k){(t>k\T j m) = 

lr ft 1 \t rC / J i • r\ 

k j=0 m=—j 


N j 

mm 2 E Pk E 2 E m Tr(T) m T) m ) |t, 


{PfcdAfc)} 


(<Pk) |2 
jm I 


k j=0 m=l 

N / j 


min 2 [ 2 m 2 pk 

{PkM} ^ ^ ^ 

j — l \ m —1 k 


\{(t>k\T 3 m\(t>k)\ 2 

ll^-mll 2 


(5.33) 


where ||2l|| 2 = (A\A). The function min renders the above equation highly 
impractical, however a useful bound can be found by recalling that pk defines a 
probability distribution which implies that 

J2Pk\(MTjm\^k)\ 2 > I ^2pk(<f>k\fjm\<l>k)\ 2 = \Tv(gf jm )\ 2 . (5.34) 

k k 

Therefore, we can drop the min function and get the final result 


N 


He) > 2 ^ 

3 =1 



= 2 1(g ), 


(5.35) 
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where (A) = Tr((5A) and the inequality is saturated only for pure states. Thus 
we have identified new quantity 1(g) (I will refer to it as quantum moment 
of inertia (QMI)) which can serve as a criterion for dynamical entanglement. 
Indeed, according to inequality (15.351) QMI is bounded for separable states: 

1 a N 

-^-(tlseparable) ^ ~cj^~( ^separable) ^ ^ 1 (5.36) 

and it reaches the maximal value for Schrodinger cat states 

m < \Hq) < ^ , 2:(|Cat v » = I.F(|Cat„» = (5.37) 

The drawback of T(q) is that it is not an efficiency criterion, since this quantity 
is not directly related to efficiency of any particular task. However, it is re¬ 
lated to dynamical susceptibility, therefore QMI is indirectly related to atomic 
interferometry. 

The QMI generalizes the analogy between dynamical entanglement and pseudo 
moment of inertia of the matrix elements of Schrodinger kitten states. Accord¬ 
ing to the Wigner-Eckhart theorem [53] the matrix element of a spherical tensor 
operator is given by 

(J, M\Tj m \J, M') = (J, M'\ j, m\J, M)^===- (5.38) 

where the reduced matrix element (J||Tj||J) is independent of m, M and M '. 
The Clebsch-Gordan coefficient (J, M'\ j, m\ J, M) does not vanish if and only if 
the condition M' + m = M is satisfied. Hence, the matrix of Tj m is sparse and 
the non-zero elements are located only in a band positioned m units away from 
the diagonal. It follows that the expectation value ( Tj m ) = Tr (gTj m ) is not zero 
only when the density matrix has a non-vanishing matrix element located in this 
band. For example, the density matrix of a kitten state (1^, M) + |-y, —M ))/\/2 
decomposes into tensors with to = ±M and to = 0 which span the diagonal. 

The quantities |(T)m)| 2 /ll^jm|| 2 appearing in the expression for 1(g) can be 
interpreted as a mass density at distance to of a token rigid body labeled by the 
quantum number j. By the analogy with the classical moment of inertia, the 
quantum moment of inertia of this body is given by the integral over distance of 
the square of to times the density (the integral turns into summation since the 
distance is quantized): ^KTjm)] 2 /\\Tj m \\ 2 ■ The total moment of the state 

is a sum of moments of all bodies - hence the summation over j in Eq. (15.351) . 

The fact that such an analogy can be established is not so surprising. The 
dynamical susceptibility inherits form the QFI its interpretation as a suscepti¬ 
bility of the state to change induced by a rotation around the z-axis (see Sec 
& The classical moment of inertia is also a susceptibility to rotations around 
given axis of the rigid body - the angular momentum of the body is proportional 
to the moment of inertia. 

Aside from this intuitive picture, the decomposition in the basis of spherical 
tensors allows to establish a relation between the dynamical entanglement of 
the state and its rotational symmetry. The state is considered more symmetric 
when the corresponding density matrix is spanned by spherical tensors with 
small angular momenta j. Indeed, if the density matrix contains only tensor 
with j = 0, then it is invariant to rotations, i.e. it is a scalar. In that case the 
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state is useless as a probe for an interferometer since no information about the 
rotation angle 9 can be imprinted. Consequently, the dynamical susceptibility 
of such a state is zero. The state is less symmetric when it decomposes into 
spherical tensors with high j-s. For example, imagine that the density matrix 
contains a maximally polarized spherical tensor Tjj (i.e. tensor with maximal 
magnetic number m = j). When the rotation around the 2 axis is applied to 
this state, the tensor transforms according to 



k times 


(5.39) 


Hence, the higher the j the more abrupt change due to rotation. Therefore, the 
less symmetric states are more susceptible to rotations which means that the 
speed along the path in density matrix space, as quantified by dynamical sus¬ 
ceptibility, can potentially reach higher values. This shows that the dynamical 
entanglement is directly related to rotational symmetry properties of the state. 


5.3 Relation between rotations and dynamical 
entanglement 

Throughout this thesis it was pointed out many times that rotations are coherent 
transformations (i.e. act on each party individually and in the same way) and 
as such cannot introduce entanglement between particles. Nevertheless, the 
example of the twin Fock state shows that rotations do affect the degree of 
entanglement detected by the dynamical susceptibility. In this section I explain 
this seemingly paradoxical behavior and provide a detailed description of the 
relation between the dynamical entanglement and rotations. This discussion 
also provides deeper insight into structure of non-classical correlations. 

In the previous section I showed that a quantum state of an fV-qubit sys¬ 
tem decomposes into parts with different symmetry properties. Each part is 
spanned by a set of spherical tensor operators, with fixed an¬ 

gular momentum j and varying magnetic numbers m. For each j I interpreted 
the expectation values of Tj m as a distribution of fictional mass over spatial 
dimension measured in units of m. This interpretation allowed for identifying 
an analogy between dynamical entanglement and the moment of inertia of a 
classical rigid body. 

The body itself can be imagined as a massive sphere which is able to rotate 
around axis passing through its center. The quantum nature of the system 
manifests itself by forcing the mass to be distributed only over certain parallels, 
such that the distance from the surface to chosen rotation axis is quantized (see 
Fig. [5711) . 
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Figure 5.1: The classical rigid body modeling the part of density matrix with 
given symmetry, spanned by spherical tensor operators with angular momentum 
j. The sphere can rotate about vertical axis passing through its center. The 
mass (red) is enforcer to only occupy parallels for which the distance from the 
axis is quantized and measured in the units of magnetic numbers to. 


The total mass around parallels \m\ units away from the axis 


equaltQ 


Total mass 

distributed over parallels 
distanced \m\ units from axis 


|( f 3m )| 2 | |( T -_ m )| 2 
\\ T jm \\ 2 || T ,_ m ||2 " 

|(f , m )| 2 + l(f^)l 2 _ J(f jm )| 2 

||2j'm|| 2 I l^jm, 11 2 


(5.40) 


Here I used the property of spherical tensor operators that Tj _ TO = (— 1 ) m Tj m . 
Therefore, the total mass of the sphere is given by the sum of masses occupying 
each allowed parallel, i.e. 


/ Total mass \ _ \ ^ / Total mass 

V of sphere ) / j \ distributed over parallel 

parallels 

Since the maximal distance form the axis to parallel at which the mass can be 
found is to = j, the radius of the sphere is given by the quantum number j. The 
quantum state is thus analogical to collection of such spheres with radii ranging 
from j = 0 to j = N (see Fig. 15.21) . 

The degree of correlations detected by the dynamical susceptibility is tied to 
the total moment of inertia of the spheres in respect to common rotation axis. 
It is important to bare in mind that according to convention adopted for the 
definition of the dynamical susceptibility, the axis of imagined rotation coincides 
with the 2 -axis of the reference frame. 


= E 


\{T 3 mW 


m=—j H-b' 1 


(5.41) 


4 Note that the magnetic number can be negative or positive, m = d=|m|. One possible con¬ 
vention is to ascribe positive/negative m to parallel located on northern/southern hemisphere. 
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Figure 5.2: The classical model of quantum state of N bosonic qubits. Each 
sphere represent the part of density matrix with distinguished rotational sym¬ 
metry (see Fig. 1,5.11) . The mass distribution one the surface of spheres (red) 
determines the moment of inertia of the system in respect to rotations about 
common axis. The moment of inertia is directly related to dynamical entangle¬ 
ment of the state. 


The correspondence between some of the traits of the classical model and 
quantum states is readily apparent. For example, an obvious feature of a clas¬ 
sical sphere is that a rotation does not change its shape, which includes total 
mass as well as radius. This is also the case for a quantum state. It was noted 
previously that rotations conserve the angular momentum j of spherical tensor 
operators, i.e. for arbitrary rotation R we have 


E ^ (K f irn&)}] = E D l',m E^’ = j(j + 1 ){RT ]m tf) ■ 

m' i 

(5.42) 


i=x,y,z 

Here 


3 

RT jm tt = E (5.43) 

m——j 


and D L'm are elements of Wigner D-matrix corresponding to R. Therefore, 
the part of the density matrix spanned by a set of tensors can 

be represented by the sphere with the same radius j even after the state has 
been transformed g RgR). The total “mass” given by Eq. (15.411) is invariant 
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as well: 


j . j (2 j + l){Tj m '){Tj m ii)* 

/ j 11 Arp £>+112 / j m' ,ra\' / “ y m" ,m) 1/iVllrn II jV\|2 






= E 


m.m' ,m" 




= E 0 '. m ' 




i(fr,iif)i 

(2j + l)<f J w)<r j -n,»)* _ y- l(^w)| 2 

' |/JV||t-i II TV \ 1 2 ll'f. 112 

KtIHHIT/I m> W-Lj ™'II 


(5.44) 


Here I used the definition of the Wigner D-matrix (see footnote[3]below Eq. (15.141) ) 
and the fact that the norm of a spherical tensor is invariant and depends only 
on j, 

\\f jm \\ 2 =Tr(T ^ m f, m ) = (2i + l)- 1 |(f||T J ||f)| 2 . (5.45) 

Hence, rotations cannot move the mass between spheres. In terms of the prop¬ 
erties of the quantum state this means that coherent operations do not change 
the symmetry of the state, as they leave unaltered the distribution of the an¬ 
gular momenta j in the decomposition of the density matrix. If the mass was 
not conserved, then it would be possible to transfer it to larger spheres which 
in turn would be equivalent to lowering the symmetry of the state. As it was 
shown previously, lower symmetry, i.e. the presence of higher j tensors in the 
decomposition of the density matrix, leads to potentially higher degree of the 
dynamical entanglement [f| 

However, it is not the total mass of spheres but rather the mass distribution 
over their surfaces that determine the degree of non-classical correlations de¬ 
tected by the dynamical susceptibility. In case of classical model, the rotations 
of spheres do not change this distribution. This can be seen from the equivalence 
of an active and a passive rotation transformation. The active transformation, 
where the point of view of an observer (represented by a choice of reference 
frame) is fixed while the rigid body is rotated, cannot be distinguished from a 
passive transformation, where the body is still and the observer is rotated. It is 
clear that the change of the reference frame has no influence on the dynamical 
properties of a body (including the mass distribution), therefore the same is true 
for rotations of the body itself. Hence the mass distribution over the surface 
of each sphere is also an invariant property of the system, independent of the 
choice of reference frame. 

The same arguments and conclusions are also valid for quantum state, with 
important modification which takes into account the quantization of magnetic 
numbers m. While the distribution of the angular momenta j in the decom¬ 
position of the density matrix depends only on the properties of the state, the 
distribution of to-s also depends on the choice of quantization axis (i.e the choice 
of reference frame). The spherical tensors Tj m are defined as the eigenopera- 
tors of the superoperator J z = [J z , .] (see Eq. (I5.23D 1. similarly to eigenstates 

5 The mass can be transferee! between spheres by, for example, non-linear evolution e~ lx L 
induced by atom-atom interactions. Such operation do not conserve angular momentum of 
spherical tensors. Therefore, interactions between particles alter the distribution of angular 
momenta in decomposition of density matrix which can decrease the symmetry of the state 
and generate additional dynamical entanglement. 
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| j, to) defined as an eigenstates of z component of the spin operator. When the 
quantization axis is changed, i.e. the operator J z is rotated into J n = RJ Z W , 
both eigenoperators and eigenstates are transformed as well: 

| j, to) -*• .R| j, to) : J n {R\j,m)) = (RJ z R^)(R\j,m)) = m(R\j,m)), (5.46) 

Tjm RTj m R) : [J n , (Rf jm R?)] = R[J z ,Tj m }R^ = m(RTj rn R^). (5.47) 

Therefore, rotation of the state, or equivalently, rotation of the reference frame 
causes displacement of allowed parallels on the corresponding classical sphere. 
This assures that the mass can be found only at quantized distances from the 
axis of imagined rotation (which by convention overlaps with the quantization 
axis), even when it shifts to new orientation. 

Although the total mass and its distribution are invariant, the dynamical 
susceptibility related to the moment of inertia calculated in respect to z-axis 
must be affected by the change of the reference frame. Indeed, as the axis 
is rotated, the distance to a given point on the surface of the sphere changes 
because the latitude of the parallel occupied by the point is now different. For 
example, consider a case when the mass is mostly distributed on the poles of 
the sphere (see Fig. 15.31 (a)). Of course, the axis passes through the poles, so 
that the moment of inertia of the sphere is very small. It is obvious that an 
appropriate rotation can reorient the z-axis so that in the new reference frame 
the mass is found on the equator instead of poles (see Fig. 15.31 fblf. The moment 
of inertia calculated for the new axis will be much larger then previously. Hence, 
even though the degree of dynamical entanglement detected for the state before 
and after the rotation will be different it is clear that the rotation itself did not 
entangle particles. Instead, it merely changed our point of view so that we could 
observe the entanglement which was already present in the state. 



Figure 5.3: According to convention adopted for dynamical susceptibility the 
moment of inertia of classical model of the state is calculate in respect to z- 
axis of the reference frame. Because of that, the dynamical susceptibility is 
affected by rotations even though the mass distribution is not. For example, in 
the case presented in (a) the mass is localized on the poles of the sphere and 
the moment of inertia is small. The same mass distribution yields significantly 
greater moment of inertia when viewed in the rotate frame (b). 

The degree of this innate dynamical entanglement existing in a given state 
g can be defined by finding the optimal orientation of the axis for which the 
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dynamical susceptibility is the largest 

( the innate \ ^ ^ ^ . 

dynamical entanglement ) = ma l xJ r (RgR J ) = J 7 (R m ax^max) (5.48) 

of state q ' R 

In case of a classical sphere the optimal axis is one of so-called principle axes 
for which the matrix of moment of inertia tensor has the largest eigenvalue. 
In this reference frame, on average, the mass is distributed closest to equator. 
Corresponding transformation of quantum state reorients quantization axis so 
that the expectation values of spherical tensors with largest magnetic numbers 
are maximized at the expense of those with smaller m- s. Note that the optimal 
axis for each sphere (and each part of the density matrix it represents) is in 
general different from the optimal axis for the whole state. The quantum state is 
represented by a collection of spheres, each having in general different orientation 
of its optimal axis. Since any rotation applied to the state affects all spheres 
simultaneously their relative orientation is conserved. Hence it is impossible to 
align the z-axis with all optimal axes of each sphere at the same time. The only 
exceptions are when the optimal axes are already aligned, or when the density 
matrix of the state is spanned exclusively by tensors with the same angular 
momentum j > 0 (j = 0 never contributes to value QMI). 

The ability to identify transformation f? max which exposes the innate dynam¬ 
ical entanglement of state is of particular usefulness for interferometry. Knowing 
the transformation allows to design interferometric sequence for a given input 
state, which results in the maximal precision of the parameter estimation. Ac¬ 
cording to Eq. (14.211) which establishes the connection between the dynamical 
entanglement and the QFI we have 

Fq(q-, Ue(R max )) = f(r max eR Lx) (5-49) 

Here the interferometric transformation is given by 

Ue = U a U n e- ieJ *UlUp, (5.50) 

and the operations U a , Up and U n which precede and follow parameter imprint 
phase satisfy condition 

-Rmax = UpUoJJn ■ (5.51) 

The formulation of the classical sphere model of the quantum state of bosonic 
qubits finalizes my inquiry of the physical interpretation of particle entanglement 
found in cold-atom systems. 


5.4 Simple example 


In this section I illustrate the ideas presented above with an example of a system 
composed of two bososnic qubits. This system is equivalent to system consisting 
one pseudo-particle with spin J = N /2 = 1. Here I will focus on the two-qubit 
variants of previously discussed states: twin Fock state, (see Eq. (15.121) 1 and the 
cat state (see Eq. (15.111) 1 


|l>L|l)« = <S[|4,i)|i,—i)] =|l,0> = |* + > 

|Cat) “ 7i " v! 


l*+> 


(5.52) 

(5.53) 
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Recall that entanglement of a state |\E' + ) is not detected by the dynamical 
susceptibility, because the z-axis is not the optimal axis (see Eq. (15.131) 1. 


,F(|T + })=4<A 2 J Z }* + 


<1,0|J 2 |1,0)-«1,0|J Z |1,0)) 5 


= 0. 


(5.54) 


In order to reveal the innate dynamical entanglement of the twin Fock, the 
reference frame has to be rotated, for example, about the a;-axis by an angle 
7t/2: 


J-( e -*2^|^ + )) = 


= 4 
= 4 


(l,0|(e^J 2 e ^ J ») 2 |l,0)-«l,0|e^^J z e *2^11,0)) 2 
(1,0|J 2 |1,0) - «l,0|J y |l,0» 2 j =4 = 7V(f + 1), 


(5.55) 


where I used the fact that e 2aJx J z e laJx = cosaJ z — sina j y . As for state |<I> + ), 
it is the maximally entangled cat states, so we have (see Eq. (14.101) 1 

J r (|4>+)) = 4 = N 2 . (5.56) 


Now I proceed to reproduce the above results with help of concepts intro¬ 
duced in the thesis and without referring to the prior experience with N > 2 
variants of the states. The first step is to decompose the density matrices of 
our states of interest in the basis of spherical tensors. In case of two-qubit, i.e. 
spin-1 system, the spherical tensors of angular momentum j = 1 are given by 
(see Eqs. (15.211) and (15.221) ): 


* - ( 1 ° ° \ 

Tio = J z = 0 0 0 , (5.57) 

V 0 0 “I / 

1 / 0 -1 0 \ 

T n = J+= o 0 -1 ) = -T\_ x . (5.58) 

v 2 \ 0 0 0 / 

The set of spherical tensor operators with j = 2 (the maximal j of basis tensors 
for J = N/2 = 1 system, see Eq. (15.291) 1 can be found by employing the Wigner- 
Eckhart theorem (15.381) and treating it as a definition of matrix elements of 
Tj m . The parameter (J = 111X211«/ = 1) sets the norm of tensor operators (see 
Eq, (15.451) 1 and can be chosen at will. By setting it to equal \J2 x j + 1 = y/5 

we get ||T 2m || 2 = 1. The explicit form of T) m matrices are given by 

1 f 1 0 °\ 

T 20 = — 0 - 2 0 , (5.59) 

v 6 \ 0 0 1 / 

1 ( 0 ^ 0 \ 

f 21 = — 0 0 1 = , (5.60) 

v 2 y 0 0 0 / 

/ o ° 1 \ 

f 22 = 0 0 0 = fJ 2 _ 2 . (5.61) 

\ 0 0 0 / 
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The decomposition in the basis is given by the standard formula for any linear 
space 


/ , 


3=0 m——j 


II T 3 


jm | 


I \t 3 


jm | 


= E E 

3=0 m=—j 


(2j m }^ + / $+ Tj, 


\\Tj 


3m 1 


IT- I 


(5.62) 

with the scalar product defined as (A|B) = Tr(A^B) (see Eq. (15.291) 1. The 
last equality comes from the direct relation between the scalar product and the 
expectation value 


(T jm \g) = Tr (f] m e) = (T jm )* s . 


(5.63) 


A straightforward calculation gives us the desired decompositions 


Q*+ = E - \/f^ 2 0 > ( 5 - 64 ) 

0 $+ = gl + ■ y ^^' 20 2 (^ 22 ~^ r ^ 2 ~ 2 ) ■ ( 5 . 65 ) 


Note that the unit matrix is spanned by j = 0 tensor, 1 oc Too- 

From Eqs. (15.641) and (15.651) we see that the density matrices of the states 
are spanned exclusively by unity and a set of tensors with angular momentum 
j = 2. Hence the states can be represented by a single sphere of radius two. 
This implies that the states have only one optimal axis, and it is possible to 
align it with the 2 -axis of the reference frame. 

The decompositions give us a general idea about the distribution of the mass 
around the spheres representing the states. The state |'F + ) is spanned only by 
m = 0 tensor, so the totality of its mass is located at the poles of the sphere. In 
case of |$ + ) state, a portion of its mass is also located the the poles. However, 
it is smaller amount then for |\f r + ), which can be seen by comparing masses of 
parallels distanced m = 0 units away from axis: 

l(T20) $+ l 2 1 l<r a o)* + | a 2 

IIT 20 II 2 6 ||f 20 || 2 3 

(recall I set the norm of j = 2 tensors to one). The rest of the mass of the sphere 
is located at the equator, which is evident from the presence of T 2 ± 2 tensors in 
the decomposition. 

We can gain a better insight into detailed distribution of mass by examining 
transformation properties of the states. For example, the quantity |(T 2 o)<p + | 2 / 
||T 20 || 2 only tells us that certain amount of mass can be found at the poles 
of sphere representing state |4/_)_). Just from this one number we cannot say 
whether this mass occupies only one of the poles or it split between them in 
some manner. However, when we examine the distribution from different point 
of view, i.e. we rotate the state, the details of the structure can become more 
apparent. In order to track the changes due to transformations more easily I 
adopt a following notation: the mass distribution on a sphere of radius j is 
represented by a vector where each row is a mass accumulated at the allowed 
parallel (with the convention that negative magnetic number refers to parallel 
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at the southern hemisphere, see footnote 0] just before Eq. (15.401) 1: 


g O 


\{ t 3] )\ 2 1 J 

I IT,, 11 2 
KT ,,-!)! 2 

!| T ,,- i || 2 


. l | T ,- j || 2 Jg 

In this notation the states 4> + and <h + are described by 


' 0 ' 

j=2 

‘ 1 " 

4 

0 


0 

2 

3 

, p$ + <-> 

1 

6 

0 


0 

0 

1 '+ 

1 

_ 4 _ 


(5.67) 


(5.68) 


Now I shall consider a sequence of two rotations about x-axis by 7 t/ 4 angle 
applied to state |4 > +). Due to transformation g<a + e~ laJx g^ + e laJa: the mass 
distribution changes according to Eq. (15.43[) which is reflected in my notation 
as 




r 1 i 

2 


0 


16 

1 


■ l “ 

4 

0 

2 

e -‘f jx 

4 

1 


0 

1 

3 


24 


6 

0 


1 


0 

0 


4 

1 


1 

. 4 . 



L 16 J 

n 



■f->- g$, 


(5.69) 


By tracing the sequence (15.691) . illustrated in Fig. 15.41 we see that initially 




Figure 5.4: Transformation of |4' + ) state as described in Eq. (MW illustrated 
on the massive sphere model: a) The initial state, b) Initial state is rotated 
through 7t/4 angle about x-axis. c) The state is rotated through additional n/4 
angle which transforms it into |•!>_}. 

the mass of state |4' + ) must have been equally distributed among the poles of 
the sphere. If this was not the case, then the first rotation could not produce 
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equally occupied northern and southern parallels distanced m = 1 units away 
from the axis. Also, since the rotation is performed about the x-axis, the mass 
is kept during the whole sequence in the z — y plane, instead of being, for 
example, “smeared” around allowed parallels. Finally, we discover that the 
second transformation, which composes with the first one to rotation by tt/2 
angle, brought state |'F + ) into |$+). Thanks to this occurrence, we also gained 
an insight into details of mass distribution of the cat state |$+): it has an equal 
fractions of its mass placed on both of the poles and what is remaining is split 
evenly on equator, at the point where the y-axis intersects the surface of the 
sphere. 

In the last step I shall utilize the decompositions (15.6411 and (15.651) of the 
states to calculate quantum moment of inertia (see Eq. (15.351) 1. which in the 
case of pure states equals the dynamical susceptibility. I start with the cat 
state | <!>+): 



(5.70) 


In case of |4/ + ) state we immediately get 

F((jisr + ) = 4X(£® + ) = 0 , 


(5.71) 


because it is spanned only by m = 0 tensors which do not contribute to QMI. 
However, we learned form examining the sequence (15.601) that state |\& + ) pos¬ 
sesses large degree of innate dynamical entanglement: 


max ?(&&»+&) = •F(e _i *- / -g# + e i *' / ") = F{e* + ) = 4 = N 2 . (5.72) 


R 


Therefore, the |'3/+) state contains the same amount of innate dynamical entan¬ 
glement as the cat state |$+). The fact that these states are equally entangled 
is unique feature of two-qubit systems. For higher numbers of particles twin 
Fock states are still significantly entangled, but to the lesser degree then the 
maximally entangled cat states. 











Chapter 6 

Conclusions 


The main purpose of the thesis was to better understand the non-classical cor¬ 
relations in ultra-cold atom systems. The main difficulty in achieving this goal 
lies in the general and discriminative definition of entanglement. One of the 
solutions to this problem, which was presented in the thesis, is to examine the 
entanglement form the point of view of its usefulness for certain tasks, like 
precise metrology. The efficiency at which this task is performed, when the 
quantum state is used as a resource, can serve as a measure of the degree of the 
entanglement. 

The atomic interferometry, a branch of quantum metrology, exploits the wave 
nature of matter and the particle entanglement to attain precise measurements 
of the phase beyond capabilities of classical devices. It turn out that ultra-cold 
atom systems are most well suited for this type of task. 

In modern experimental realizations of such interferometers the Bose-Einstein 
condensate (BEC) used as a probe can be described as a collection of qubits. 
Due to their bosonic nature, the state of the BEC can be described as a state 
of a single pseudo-particle with the spin equal to the half of its particle number 
N. It follows that any (unitary) interferometric transformation is mapped onto 
a sequence of rotations of the spin, with one of the angles being the unknown 
phase. Moreover, it is important to note that all possible unitary transforma¬ 
tions of a single qubit are equivalent to rotation, and it follows that any coherent 
transformation of the collection of qubits (i.e. the transformation which acts on 
each qubit in the same way) is also mapped onto global rotation of the system. 
Since the indistinguishable qubits cannot be addressed individually, coherent 
rotations constitute the entirety of possible local operations - transformations 
which do not introduce or destroy the entanglement between particles. 

The efficiency of the interferometer is tied to the precision of the phase esti¬ 
mation. According to the Cramer-Rao theorem, which is known from classical 
theory and was later supplemented by Braunstein and Caves with quantum con¬ 
siderations, the precision is bounded by the quantum Fisher information (QFI). 
The QFI quantifies the susceptibility of the state to change induced by a given 
transformation, which in this case is the interferometric sequence. 

When the transformation is coherent, the QFI cannot surpass the threshold 
of the shot-noise limit unless the state is entangled. Therefore, for a given 
interferometric sequence, the QFI can serve as a criterion for the entanglement 
that is useful for enhancing the precision of this particular interferometer. For 
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another choice of the interferometer, the QFI would deem different set of states 
as usefully entangled. 

On one hand, the freedom of choice of the transformation is essential for 
embedding the question of entanglement nature in the context of a real experi¬ 
ment. On the other hand, this ambiguity might obscure some of the important 
structures and relations which could help in understanding the non-classical 
correlations, which is the main aim of this thesis. This problem was solved 
by introducing the dynamical susceptibility, which is defined as the QFI for 
“standard” interferometer and is a function of the state alone. The dynami¬ 
cal susceptibility allowed to show that states which are usefully entangled for 
one interferometric sequence are related to states useful for different interfer¬ 
ometer through coherent transformation, which conserves the degree of particle 
entanglement. 

By exploiting the equivalence between bosonic qubits and the spin system, it 
has been found how particle indistinguishability can enhance the overall degree 
of the entanglement. The susceptibility to change due to coherent transforma¬ 
tion (quantified by the dynamical susceptibility) of the spin N/2 system is in 
general greater than the susceptibility of the ensemble of N individuals with 
spins 1/2. If the qubits constituting the state were distinguishable, the spins of 
the individuals could add-up to a whole range of total angular momenta, thus 
diminishing the state’s susceptibility. For identical qubits, the only possibility 
is to add the individual spins to the maximal total spin of N/2 - allowing for 
potentially maximal degree of entanglement. 

The entanglement criteria based on the QFI, including the dynamical sus¬ 
ceptibility, establish a direct correspondence between the particle entanglement 
and the susceptibility of the state to rotations. In classical physics such suscep¬ 
tibility is measured by the moment of inertia of the system, which is given by 
the distribution of the mass in respect to the axis of rotation. It was shown that 
the dynamical susceptibility is indeed closely related to the quantum analog of 
the classical moment of inertia. In the quantum case the axis of rotation is given 
by the quantization axis and the classical mass is replaced by the components 
of the density matrix of the state decomposed in the basis of spherical tensor 
operators. 

The final result of the thesis is the formulation of the classical model of the 
quantum state and its entanglement. Within this model the state of N bosonic 
qubits is represented by nested massive spheres of increasing radii. The quantum 
nature of the system is manifested by enforcing the mass of each sphere to 
only occupy parallels located at quantized distances from the common rotation 
axis. The model embraces the analogy between the dynamical susceptibility 
and the classical moment of inertia of a rigid body and identifies the degree of 
correlations between qubits with the distribution of the fictional mass on the 
surfaces of the spheres. The most entangled states are represented by the system 
with the majority of its mass located at the equator of the largest sphere so that 
the moment of inertia is the greatest. 

I hope that the existence of the model such as this can convince the Reader 
that at least certain type of particle entanglement can be understood in terms 
of concrete physical properties of the system. That the difference between the 
classically-correlated state and the entangled one is only quantitative not qual¬ 
itative. There is no sharp border between the two cases: separable states are 
characterized by a lesser “intensity” of some special physical property then the 
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states that are considered as entangled. I believe that it is only natural to ex¬ 
pect such smeared border rather then the sharp one. After all, in the quantum 
world there are no “classical correlations” because everything is just “quantum”. 
Hence, the distinction between separable and entangled state is in most cases 
artificial, and so no spectacular physically effects should accompany the passing 
of the border. 
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